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the  characteristics  of  the  observed  fluctuations  are  consistent  with 
the  theory  of  linear,  equatorially  trapped  waves,  particularly  the 
Yanai  wave. 

Numerical  integration  of  equations  of  this  system  shows  that 
oscillations  such  as  those  observed  may  be  a resonant  response  to  an 
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1 . INTRODUCTION 


The  equatorial  regions  of  the  oceans  and  atmosphere  have  long  been 
a topic  of  intense  interest  to  the  theoretician  as  well  as  to  the  obser- 
vationalist.  Unique  aspects  of  the  behavior  of  the  circulation  in  these 
regions  are  primarily  the  result  of  the  zero-crossing  of  the  vertical 
component  of  the  Coriolis  acceleration. 

At  latitudes  away  from  the  equator,  circulation  problems  can  be 
treated  with  quasi-geostrophic  theory.  This  theory  is  a manifestation 
of  the  constraints  on  long-period  motions  resulting  from  the  earth's 
rotation.  The  gross  aspects  of  the  circulation  can  be  understood  with- 
out inclusion  of  non-linear  and  complicated,  frictional  dependence. 

Near  the  equator,  this  is  not  the  case.  It  is  not  immediately 
apparent  which  balances  are  important.  Quasi-geostrophic  theory  takes 
on  a different  interpretation  since  the  inertial  period  is  very  mu'h 
dependent  on  the  latitude  and  becomes  infinite  at  the  equator.  While 
the  theory  can  be  defined  using  an  equatorial  adjustment  time  scale 
( Moore  and  Philander,  1976  ),  there  is,  nevertheless,  a latitude 
at  which  the  Coriolis  acceleration  becomes  so  small  as  to  result  in 
balances  which  include  non-linear  and  dissipative  processes  in  the  steady 
state  . 

On  the  other  hand,  time-dependent  motion  could  exist  in  which  non- 
linear and  dissipative  processes  are  not  significant.  This  would  be 
true  mainly  in  the  absence  of  currents  with  intense  shear,  which  is 
not  the  case  in  the  real  ocean. 

The  equatorial  ocean  is  characterized  by  intense  zonal  currents 
to  a depth  of  about  400  m.  These  currents  display  strong  variations 
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in  space  and,  based  on  more  recent  observations,  time.  Typically , from 
the  surface  to  a depth  of  about  50  m,  from  about  5°N  to  5°S  in  the 
Atlantic,  there  is  westward  flow  with  speeds  up  to  120  cm/sec.  From 
about  50  m to  200  m or  more,  the  flow  is  eastward  with  a somewhat  smaller 
magnitude.  Below  this  depth,  there  is  sometimes  a weaker  westward  flow. 
The  intense  eastward  flow  is  referred  to  as  the  Equatorial  Undercurrent 
or  Lomonosov  Current.  In  the  Pacific,  a similar  structure  exists 
(Philander , 1 973b ; Knauss,  1966;  Stalcup  and  Metcalf,  1966;  Rinkel,  1964; 
and  Knauss,  1961). 

Variations  from  periods  of  several  months  to  a year  have  been  noted 
in  earlier  studies,  although  evidence  of  h igher- frequency  fluctuations 
has  been  presented  by  Rinkel  (1964),  Wunsch  and  Gill  (1976),  and  Groves 
and  Miyata  (1967). 

During  the  GARP  Atlantic  Tropical  Experiment  (GATE),  overwhelming 
evidence  was  found  for  the  existence  of  energetic  fluctuations  of  equa- 
torial fields  with  periods  ranging  from  tidal  to  several  months  (Duing, 
1975).  The  result  has  been  an  increased  interest  in  time-dependent, 
equatorial  motion.  There  was  also  evidence  of  energetic  fluctuations 
in  the  tradewinds  at  similar  periods  (Kr ishnamur t i , 1975  a and  b). 

It  is  not  clear  from  the  observations  alone  whether  the  current  fluctua- 
tions constitute  a true  meander  or  whether  they  are  part  of  a wavelike 
process  which  exists  over  a wider  latitude  range. 

According  to  Moore  and  Philander  (1975)  and  Gill  (1975),  theories 
on  equatorial  motions  can  be  divided  into  two  somewhat  exclusive  cate- 
gories: linear  and  non-linear,  steady-state,  wind-driven  models  (mostly 

numerical),  and  linear,  inviscid,  time-dependent,  free  and  forced  models 
(mostly  analytical).  Notable  exceptions  include  studies  by  Cane  (1975) 
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and  by  O'Brien  and  Hurlburt  (1973),  where  time  dependence,  advection, 
and  friction  have  been  included.  Even  these  latter  investigations, 
however,  are  primarily  concerned  with  the  long-term  spinup  of  an  entire 
ocean  initially  at  rest,  rather  than  short-term  variability. 

Non-linear,  steady-state  models f which  are  typified 
in  a paper  by  McKee  (1973),  are  in  essential  agreement  as  to  the  reasons 
for  the  existence  of  the  Equatorial  Undercurrent  and  the  general  equa- 
torial circulation.  They  predict  correct  scales  in  latitude  and  zonal 
velocity  and  agree  on  the  gross  features  of  the  meridional  or  secondary 
circulation.  They  rely,  to  a large  extent,  on  vertical  and/or  horizontal 
eddy  viscosity  coefficients,  however,  and  thus,  finer  structure  in  their 
results  is  somewhat  arbitrary.  The  vertically  averaged,  non-linear 
terms,  in  conjunction  with  the  meridional  circulation,  tend  to  transport 
eastward  momentum  equatorward  (Gill,  1975).  A common  interpretation 
of  all  steady-state  models  is  that  an  easterly  wind  over  an  equatorial 
ocean  bounded  to  the  west  will,  within  several  weeks  after  onset,  develop 
an  eastward  pressure  gradient  force  as  a result  of  set-up  against  the 
boundary.  The  setup  is  the  result  of  westward  transport,  which  occurs 
primarily  in  the  surface  layer  near  the  equator.  An  eastward  flow  results 
at  a depth  where  this  pressure  gradient  force  exceeds  the  westward  drag 
of  the  wind,  A secondary  circulation  develops  due  to  the  poleward  Ekman 
divergence  in  the  mixed  layer,  resulting  in  an  equatorward  flow  in  the 
thermocline.  This  meridional  flow  modifies  the  zonal  velocity  profile 
through  advective  effects. 

A difficulty  with  the  homogeneous  models  is  that  a vigorous  upwel- 
ling  develops  at  the  equator  as  part  of  the  meridional  circulation. 

This  is  not  consistent  with  observations  according  to  Moore  and  Philander 
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(1975).  In  the  presence  of  stratification,  however,  a depth-dependent 
zonal  flow  can  approach  geostrophy,  and  thus,  the  meridional  circulation 
will  be  less  than  in  the  homogeneous  case.  Stratified,  non-linear, 
steady-state  models  (Robinson,  1966;  Philander,  1973)  show  the  same 
essential  features  as  the  homogeneous  models  near  the  surface,  except 
for  the  reduced  upwelling.  Due  to  the  additional  thermohaline  effects, 
however,  the  stratified  models  help  to  explain  the  deep  westward  flow 
observed  by  Stalcup  and  Metcalf  (1966). 

While  the  steady-state  theories  cannot  explain  the  dynamics  of 
the  observed  fluctuations,  they  set  the  stage  for  possible  explanations 
which  are  offered  in  this  study.  Theory  provides  a mechanism 
by  which  the  equatorial  current  system  can  come  into  existence  as  a 
nearly  geostrophic  zonal  flow.  The  meridional  flow  constitutes  a second- 
ary circulation  which  is  a result  of  the  Ekman  divergence  in  the  surface 
layer.  In  the  steady-state,  momentum  and  energy  input  by  the  wind  is 
balanced  by  frictional  dissipation  through  vertical  and  horizontal  shear. 

Mikhailova  and  Shapiro  (1972)  have  treated  a problem  similar  to 
that  of  Charney  and  Spiegel  (1971),  except  that  a two-layer  density 
field  is  included.  The  model  also  includes  a linear  friction  law  (de- 
pending on  horizontal  velocity)  in  lieu  of  horizontal  viscosity.  They 
approach  the  time-dependent  aspect  by  applying  a time-varying  zonal 
wind,  which  results  in  an  increase  of  the  mean  transport  due  to  the 
non-linearities.  The  v-component  of  the  response  is  antisymmetric 
about  the  equator,  however,  and  the  shortest  period  in  rhe  forcing 
is  two  months,  thus  limiting  the  applicability  of  the  model  to 
fluctuations  observed  during  GATE. 
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f Philander  (1976)  has  found  that  surface  current  systems  such  as 

the  Atlantic  South  Equatorial  Current  and  the  North  Equatorial  Counter- 
current  may  be  unstable  to  certain  waves,  and  that  the  Undercurrent 
may  respond  passively  through  vertical  coupling.  The  parameters  in 
the  theory  include  the  characteristics  of  the  mean  westward  flow  on 
a latitude  scale  of  about  10°,  which  is  beyond  the  range  of  the  PCM 
GATE  data.  This  type  of  theory,  which  is  a linearized  wave  problem 
(i.e.,  with  the  mean  flow  included  as  a space-variable  coefficient), 
states  that  infinitessimal  disturbances  to  a mean  state,  provided  they 
have  the  right  characteristics,  will  grow  in  time,  deriving  energy  from 
the  mean  flow  through  non-linear  interaction.  Once  the  perturbation 
has  grown  sufficiently,  the  dynamics  must  move  from  a linearized  to 
a truly  non-linear  regime.  On  the  other  hand,  with  a finite  perturba- 
tion such  as  a time-dependent  wind  forcing  with  characteristics  that 
do  not  necessarily  satisfy  the  instability  conditions,  a response  invol- 
ving non-linear  interactions  would  also  result  which  might  be  related 
to  the  meanders  observed. 

Cane  (1975),  in  a non-linear,  time-dependent,  numerical  model*, 
includes  a case  where  a southerly  wind-component  is  added  to  the  forcing 
after  the  model  is  integrated  to  an  equilibrium  with  an  easterly  wind. 


*Cane's  model  is  composed  of  an  active  layer  of  one  density  overlying 
another  of  higher  density.  The  lower  layer  is  quiescent,  or  motionless. 
The  active  layer  is  divided  into  a surface  frictional  layer  of  fixed 
thickness  and  a thicker  subsurface  layer  which  contains  the  Undercurrent. 
Vertical  advection  of  momentum  between  the  layers  is  included,  as  well 
as  vertical  and  horizontal  viscosity. 
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A result  is  the  upwind  shift  of  the  Undercurrent  over  a period  of  about 
two  months  and  the  development  of  a quasi-stationary  meander  pattern  in 
latitude  and  longitude,  originating  at  the  western  boundary  and  propa- 
gating (in  the  group  sense)  eastward.  During  the  first  month,  over  the 
whole  domain,  especially  in  the  eastern  quarter  where  longitudinal 
structureis  initially  negligible,  a transient,  north-south  oscillation 
of  the  Undercurrent  occurs  with  a period  of  about  16  days.  It  even- 
tually damps  out,  but  this  may  be  due  to  the  dissipation  in  the  model. 
It  appears  to  be  a locallv  excited  transient.  A deficiency  of  this 
calculation,  as  far  as  CATE  observations  are  concerned,  is  that  the 
surface  zonal  flow  which  exists  at  the  onset  of  the  southerly  wind  is 
eastward  instead  of  westward  over  most  of  the  domain  along  the  equator. 
This  is  likely  related  to  the  peculiar  nature  of  the  surface  layer  in 
the  model.  These  results,  nevertheless,  bear  the  closest  resemblance 
to  the  observations  insofar  as  the  meanders  are  concerned. 

The  problem  of  linear  equatorial  waves  has  been  dealt  with  in 
detail  by  Matsuno  (1966),  Moore  (1968),  Lighthill  (1969),  and  Cane  and 
Sarachik  (1976).  The  most  extensive  study  for  free  equatorial  waves 
in  the  ocean  is  given  by  Moore  (1968). 

The  scales  of  motion  predicted  by  these  theories  are  consistent 

with  the  features  of  GATE  observations,  but  the  assumption  of  linearity 

is  not  completely  justified.  Moore  (1968)  defines  an  equatorial  Rossby 
2 

number,  Rq  = u/gj*’.  where  U is  the  velocity  scale,  L is  a length  scale, 
and  8 is  the  latitudinal  derivative  of  the  Coriolis  parameter.  It 
can  be  said  that  this  is  the  Rossby  number  at  latitude  L for  length 


scale  L. 
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If  L is  sufficiently  large  or  U sufficiently  small,  R^  may  be 
neglected,  resulting  in  the  linear  equations  of  motion.  It  turns  out 
a posteriori  that  L is  not  arbitrary,  but  constitutes  a trapping  scale 
for  equatorial  waves.  According  to  Moore  (1968), 


L = 


where  is  a gravity  wave  speed  associated  with  the  vertical  mode  under 
consideration  and  is  about  200  cm/ sec  for  low-order  baroclinic  modes. 

L,  which  is  also  termed  the  equatorial  radius  of  deformation,  R^  , is 
about  100  to  300  km  in  these  cases.  Consequently,  for  U > 50  cm/ sec, 

R^  is  no  longer  negligible.  In  particular,  if  a mean  zonal  flow,  such 
as  that  observed  near  the  equator,  is  included,  the  Rossby  number  exceeds 
unity,  regardless  of  the  velocity  scale  of  the  fluctuations. 

Linear  wave  theory,  therefore,  would  be  applicable  for  small  ampli- 
tudes if  there  were  no  equatorial  currents.  In  the  real  ocean,  however, 
these  currents  exist,  and  some  manifestation  of  their  effect  must  be 
included  in  the  study  of  time-dependent  motions  near  the  equator. 

While  the  steady-state  and  large-scale,  time-integrated  numerical 
models  include  features  such  as  the  Undercurrent,  it  is  the  wave  theories 
which  are  most  directly  applicable  to  motions  such  as  those  observed 
during  GATE.  The  question  still  remains  as  to  the  origin  of  the  distur- 
bances. 

The  instability  mechanism  proposed  by  Philander  (1976)  offers  one 
explanation.  However,  Philander  (personal  communication)  has  stated 
that  several-week  oscillations  with  wavelengths  larger  than  1500  km 
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are  not  likely  to  be  unstable  in  the  Atlantic,  and  the  observed  fluc- 
tuations appear  to  be  of  the  order  of  2500  km  in  length  (Dlling  et^  al , 
1975). 

Waves  generated  in  a different  location,  e.g.,  near  a boundary, 
which  propagate  into  the  area  of  observation,  provide  another  possibi- 
lity. Problems  of  this  nature  have  been  treated  by  Moore  (1968)  in 
the  absence  of  a mean  zonal  current. 

A third  explanation,  which  is  not  entirely  exclusive  of  the  second, 
is  wind  forcing.  This  problem  has  been  treated  formally  by  Lighthill 
(1969)  and  Cane  and  Sarachik  (1976)  for  the  equatorial  ocean.  They 
suggest  that  many  equatorial  modes  can  be  excited,  depending  on  the 
form  of  the  forcing  process.  These  studies  are  quite  general  in  scope 
and  do  not  treat  the  specific  cases  of  interest  in  any  detail.  More- 
over, they  do  not  include  consideration  of  mean  currents. 

Observational  evidence  for  wind-forced  equatorial  waves  is  offered 
by  Wunsch  and  Gill  (1976)  from  sea- level  fluctuations  in  the  equatorial 
Pacific.  They  found  that  correlation  between  atmospheric  pressure  and 
sea-level  indicated  the  excitation  of  equatorially  trapped  inertia- 
gravity  waves  in  the  first  baroclinic  mode. 

The  above  description  of  Cane's  1975  model  also  suggests  that 
several-week  oscillations  may  be  transients  excited  by  meridional,  time- 
dependent  wind  forcing.  The  advantage  of  this  model  is  that  an  under- 
current exists;  i.e.,  non-linear  effects  have  been  included. 

The  purpose  of  the  present  study  is  to  try  to  determine  the  nature 
of  the  current  fluctuations  observed  during  GATE  on  the  basis  of  their 
kinematic  and  dynamic  properties,  and  to  provide  a plausible  hypothesis 
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for  their  generation.  In  particular,  it  is  shown  that  such  fluctuations 
are  likely  to  be  a resonant  response  to  a broad-band,  impulsive  exci- 
tation rather  than  a direct  passive  response  to  a continuing  forcing 
component  or  the  result  of  a dynamic  instability  process  in  which  wave 
characteristics  depend  on  details  of  mean  flow  profiles.  It  is  further 
demonstrated  that  a meridional  windstress  component,  such  as  that  observed 
during  GATE,  is  adequate  to  produce  such  a response. 
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2.  ANALYSIS  OF  THE  GATE  EQUATORIAL  PROFILING  EXPERIMENT  (GEPE) 

2 . 1 Overview 

During  the  GARP  Atlantic  Tropical  Exper iment (GATE)  in  1974,  a great 
number  of  observations  were  obtained  in  the  Equatorial  Atlantic  by  many 
different  institutions.  It  will  be  some  time  before  all  data  are  pro- 
perly assimilated  and  combined  to  give  an  overall  picture  of  the  pro- 
cesses in  the  atmosphere  and  ocean  during  the  period  of  the  experiment. 
Descriptions  of  the  oceanographic  measurements  appear  in  GATE  Report 
#14  (1976). 

The  present  investigation  is  based  on  the  analysis  of  GEPE  data, 
surface  wind  data  from  ships  in  the  vicinity  of  GEPE,  the  results  of 
the  reports  indicated  above,  and  several  publications  by  GATE  inves- 
t igator s . 

The  GEPE  data  were  acquired  by  workers  from  the  University  of  Miami 

**  V»  * ' 

from  four  ships  during  Phase  II  of  GATE,  from  7 July  to  19  August. 

On  R/V  COLUMBUS  ISELIN,  Dlling  (1975b)  acquired  65  profiling  current 
meter  (PCM)  stations,  44  STD  stations,  and  166  XBT  traces  while  tra- 
versing a trapezoidal  track  between  1.75°N  and  1.25°S,  1 8°W  and  29°W. 
Brown  (1975)  acquired  30  PCM's  from  the  GDR  vessel  ALEXANDER  VON  HUMBOLDT 
between  1 . 5°N  ana  1.5°S,  23.5°W  and  25°W.  Hallock  (1975)  obtained  43 
PCM's  from  the  U.S.S.R.  vessel  AKADEMIK  KURCHATOV  at  23.5°W.  At  10°W, 

13  PCM's  were  obtained  by  Low  (1975)  from  the  French  ship  CAPRICORNE. 

My  analysis  efforts  centered  on  the  PCM  data  from  R/V  ISELIN  since 
they  are  the  most  comprehensive  data  from  GEPE.  Similar,  but  less  exten- 
sive, analysis  was  performed  on  PCM  data  from  R/V  HUMBOLDT.  Data  from 
R/V  KURCHATOV  and  R/V  CAPRICORNE,  being  more  restrictive  in  scope,  have 
been  used  for  comparison. 
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Surface  wind  observations  have  been  investigated  to  estimate  the 
spatial  and  temporal  variability  of  atmospheric  forcing  on  the  ocean 
in  the  GEPE  area . 

The  analysis  was  performed  in  five  main  steps:  1)  PCM  data  were 

calibrated,  edited  and  corrected  for  systematic  errors  due  to  observa- 
tional procedures;  2)  interpolation  and  smoothing  in  latitude  and  time 
were  effected  through  a hand  analysis  followed  manual  digitization; 

3)  diagnostic  calculations  on  the  analyzed  data  were  performed  to  esti- 
mate magnitudes  of  accelerations  and  momentum  balances;  4)  correlations 
between  the  surface  wind  and  current  were  computed;  and  5)  the  data 
were  averaged  into  two  active  layers  consistent  with  principal  vertical 
distr ibutions . 

The  objectives  of  the  analysis  described  in  the  remainder  of  this 
chapter  are  to  describe  the  major  features  of  the  PCM  data  set,  obtain 
quantitative  estimates  of  the  dynamic  processes  implied,  and  to  state 
a hypothesis  for  explaining  the  observed  fluctuations  of  the  equatorial 
currents . 

2 . 2 PCM  Data  Analysis 

2.2.1  Data  Acquisition  and  Analysis 

The  PCM  data  from  R/V  ISEblN  are  now  examined.  Figure  1 is  a plot 
of  PCM  station  positions  distributed  in  latitude  and  time  (y-t  plane). 
This  figure  represents  the  sampling  pattern  for  the  profiles  where 
profiles  from  29°W  are  lumped  together  with  those  from  28°W,  assuming 
that  east-west  scales  are  much  greater  than  1°  longitude. 
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At  each  position,  a PCM  and  an  STD  profile  were  made  to  a depth 
of  about  600  m.  In  a similar  manner,  XBT  profiles  were  obtained  over 
the  same  sampling  domain.  Preliminary  processing  of  the  PCM  data  was 
performed  aboard  ship  on  a PDP-8E  mini-computer,  allowing  for  early 
analysis  of  gross  features  (DUing  £t  a_l.  , 1975).  Subsequent  processing 
and  analysis  was  done  on  a UNIVAC  1106  in  Miami.  Details  of  the  data 
reduction  can  be  found  ir.  Brown  et^  al . ( 1977  ) . 

Due  to  irregular  spacing  and  paucity  of  observations  in 
y and  t,  I decided  to  do  a hand  analysis  (not  unlike  classical  meteoro- 
logy techniques)  of  the  velocity  components  for  each  10-meter  level 
to  a depth  of  200  m.  Below  200  m,  currents  become  weak  and  variable, 
and  the  results  of  this  type  of  analysis  become  questionable.  A similar 
analysis  was  done  for  the  surface  wind  and  thermocline  depth,  as  observed 
from  R/V  ISELIN.  The  resulting  contour  maps  are  typified  in  Figs.  2 to  5. 
The  maps  were  first  sketched  with  the  positions  of  the  contours  being 
estimated  by  manual  linear  interpolation  between  data  points  in  the 
y-t  plane.  Details  on  scales  less  than  the  average  y-t  spa  ing  of  the 
data  were  smoothed  out  and  derivatives  (horizontal)  were  made  reasonably 
uniform  in  areas  of  large  magnitudes.  The  maps  were  subsequently  digi- 
tized to  a one-day,  15-mile  grid  for  further  analysis.  A series  of 
meridional  (latitude-depth)  sections  was  prepared  from  the  digitized 
data.  These  sections  (Figs.  6 to  10)  are  contour  plots  of  zonal  velocity 
and  vector  plots  of  meridional  and  vertical  velocity  (more  will  be  said 
later  about  the  calculation  of  vertical  velocity).  Figures  17  a and  b 
are  y-z  sections  of  time-averaged  zonal  and  meridional  velocity  compo- 
nents (u,  v).  These  figures  show  that  the  structure  of  the  flow  field 
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is  quite  complex  and  that  it  changes  significantly  with  time,  but  that 

in  the  mean,  v is  quite  small.  It  can  also  be  seen  from  Fig.  17b  that 

a calculation  of  mean  verti-al  velocity  from  v would  not  be  signifi- 

y 

cantly  different  frora  zero  in  any  part  of  the  domain  (y,  z)  of  the 
observat ions . 

Two  properties  of  the  GATE  PCM  data  set  indicate  the  use  of  caution 
when  interpreting  results  of  time  averages.  These  are  the  shortness 
of  the  record  compared  to  the  period  of  the  principal  fluctuations  and 
the  amplitude  of  these  fluctuations  compared  to  the  mean.  In  addition, 
there  is  evidence  of  significant  fluctuations  with  periods  longer  than 
the  record  length:  in  our  observations,  the  mean  position  of  the  Under- 
current is  north  of  the  equator,  in  contrast  to  results  at  23.5°W  (Bubnov, 
1975),  as  well  as  theoretical  predictions  which  place  the  Undercurrent 
upwind  of  the  equator  in  the  mean  (eg.,  McKee,  1973;  Charney  and  Spiegel, 
1971;  Philander,  1973).  In  short,  while  we  can  get  good  estimates  of 
the  scales  of  the  time-dependent  behavior  of  the  equatorial  currents, 
mean  values  cannot  be  expected  to  be  representative  of  long-term  mean 
conditions. 

2.2.2  Description  of  Analyzed  Data. 

The  most  outstanding  feature  in  Fig.  3 is  a lateral,  latitudinally 
antisymmetric  fluctuation  in  time  of  the  zonal  velocity  component  at 
a depth  of  80  m.  Examination  of  the  meridional  component  reveals  a 
north-south  advection  consistent  with  this  fluctuation  (this  point  will 
be  systematically  demonstrated  in  Section  2.2).  The  westward  flow  at 
the  surface,  while  not  as  well  defined,  exhibits  variability  on  the 
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same  scales  (Fig.  2).  The  possibility  that  the  observed  2 to  3-week 
oscillation  is  an  artifact  of  the  sampling  scheme  has  been  ruled  out 
since  it  is  present  in  several  independent  sets  of  observations. 

The  meridional  sections  (Figs.  6 to  10)  exhibit  a distinct,  bi- 
modai  structure:  the  north-south  oscillation  of  the  westward  surface 

flow  maximum  tends  to  be  about  90°  out  of  phase  with  that  of  the  core 
of  the  Undercurrent.  This  feature  is  suggestive  of  the  mixed-mode 
oscillation  of  a two-layer  system. 

The  y-t  maps  of  the  zonal  component  at  80  m (Fig.  3)  show  maxima 
symmetric  in  latitude  about  the  position  of  the  core  of  the  Undercurrent. 
These  may  be  associated  with  symmetric  wave  modes  or  may  be  the  eastward 
acceleration  which  would  result  from  a poleward  advection  of  the  equa- 
torial water. 

Figures  11  and  12  present  the  corresponding  y-t  maps  at  10  m and 

at  50  m for  the  zonal  and  meridional  components  of  velocity,  as  observed 

from  R/V  HUMBOLDT  at  23.5°W.  By  comparing  mer id iona 1 ( v ) components 

at  10  m from  23.5°w  and  28°W,  we  can  infer  a westward  phase  propagation 
of  about  150  cm/sec,  assuming  that  the  phases  in  the  two  plots  are  one 

and  the  same.  The  calculation  is  done  by  dividing  the  difference  in 

longitude  by  the  temporal  offset,  which  is  about  4.5  days.  The  resulting 

speed  is  in  approximate  agreement  with  a similar  calculation  by  Dlling 

et  al . (1975),  who  compared  y-t  maps  of  maximum  zonal  (u)  components 

from  the  two  ships.  They  obtained  a time  difference  of  three  days, 

giving  a phase  speed  of  190  cm/sec.  The  u-component  at  the  surface 

is  not  well  enough  defined  to  identify  phases  with  any  certainty,  and 

the  v-components  at  the  core  depth  are  difficult  to  compare  due  to  a 

non-zero  mean  (which  may  be  in  error  - see  Fig.  21)  at  the  HUMBOLDT 
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longitude.  DUing  et_  aK  ( 1975)  also  estimated  a period  of  16  days  by 
identifying  quasi-sinusoidal  appearance  of  the  map  of  u-max  with  one 
cycle  of  a fluctuation  or  "meander." 

Clearly,  it  is  difficult  to  define  a period  or  a phase  speed  with 
great  certainty.  We  should,  therefore,  infer  a range  of  such  parameters 
from  the  data.  According  to  DUing  (1975),  the  best  estimate  of  the 
time  lag  between  phases  at  28°W  and  23.5°W,  based  on  salinity  observa- 
tions from  R/V  HUMBOLDT,  R/V  I SELIN,  and  R/V  TRIDENT,  as  well  as  PCM 
data,  is  about  three  days.  The  estimates  made  with  the  v-components 
at  the  surface,  as  described  above,  yield  a value  of  4.5  days.  The  east- 
west  distance  between  ISELIN  and  HUMBOLDT  is  500  km.  We  thus  obtain 
a phase  speed  ranging  from  130  to  180  cm/sec . In  all ‘cases  where  a 
qua  si-periodic  structure  in  y-t  can  be  defined,  the  estimate  of  the 
period  varies  between  15  and  17  days.  We  can  now  derive  a range  of 
values  for  an  east-west  wavelength. 

C =s  (130  to  190)  cm/sec 
T as  (15  to  17)  days 
\ « (2070  to  2790)  km 

It  must  be  emphasized  that  the  independent  observed  parameters 
are  C and  T,  whereas  \ is  derived  from  them,  assuming  a dependence 
of  the  form 

i (ot  + kx)  where  k=  2tt/\ 

Alternative  interpretations  are  possible  but  less  likely.  Following 
DUing  et  al_.(1975),  we  construct  a table  (c.f.  Table  1)  of  C,  \ , vs.  n 
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for  eastward  and  westward-propagating  waves.  The  number  of  complete 
cycles  of  the  waves  that  exist  between  the  two  longitudes  is  represented 
by  n.  Previous  calculations  have  been  for  the  n = 0 case. 


Table  1.  Wave  Parameters  Estimated  from  Observations 


n 

Cn  (cm 

/ sec  ) 

n ( km) 

i 

4!  - 

55 

607  - 714 

for  eastward 

2 

19  - 

23 

274  - 294 

waves 

0 

130  - 

190 

2070  - 2790 

1 

27  - 

32 

397  - 415 

for  westward 
waves 

2 

15  - 

17.5 

222  - 227 

This  table  summarizes  the  possible  wave  parameter  ranges  which 
can  be  derived  by  visual  inspection  from  the  y-t  maps.  These  ranges 
will  be  used  in  Chapter  4,  where  comparison  with  theoretical  dispersion 
curves  is  made. 

In  addition  to  the  y-t  maps  of  current  velocity  components,  a 

similar  map  of  the  depth  of  the  21°C  isotherm  ( thermocline)  is  presented 

in  Fig.  5.  A quasi-per iodic  structure  appears  in  this  quantity  as  well. 

The  variations  toward  the  end  of  the  record,  if  interpreted  as  vertical 

-2 

motion,  imply  oscillatory  vertical  velocities  of  theorder  of  10  cm/sec. 
This  value  is  in  agreement  with  estimates  of  vertical  velocity  computed 
from  the  divergence  of  horizontal  velocity.  Figure  5 shows  the  most 
pronounced  change  in  height,  from  55  m to  95  m,  between  10  August  and 
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14  August,  indicative  of  an  8-day  period  if  interpreted  as  a half-cycle. 
This  result  is  in  agreement  with  those  derived  from  temperature  obser- 
vations by  Miller  ( 1977)  and  Bubnov  et^  al_.  ( 1975). 

While  there  exist  many  interesting  features  in  the  observations  as 
presented,  e.g.,  small-scale  vertical  structure  fluctuations  which  might 
be  indicative  of  the  presence  of  internal  waves,  the  focus  of  this 
investigation  is  on  the  north-south  oscillations  of  the  eastward  Under- 
current and  the  westward  surface  flow.  For  purposes  of  simplification, 
the  digitized  y-t  map  data  were  averaged  into  two  layers  between  0 and 
200  m The  boundary  between  the  layers  was  chosen  as  the  depth  of  zero 
zonal  velocity  component.  One  reason  for  this  choice  is  that  it  is 
reasonable  to  expect  the  dynamics  of  a westward  flow  to  be  significantly 
different  from  that  of  an  eastward  flow  at  the  equator.  Another  reason 
is  that  the  layered  configuration  chosen  allows  a more  direct  comparison 
of  the  observations  with  a two-layer  dynamical  model,  such  as  those 
previously  discussed,  and  an  additional  one  to  be  introduced  in  Chapter 
3. 

The  y-t  maps  of  these  layer-averaged  data  are  presented  in  Figs. 

13  to  16.  It  is  seen  that  the  essential  features  of  the  lateral  oscil- 
lations are  retained,  while  minor  features  are  smoothed  out  and  ampli- 
tudes are  decreased. 

2 . 3 Dynamics  and  Kinematics 

2.3.1  Accelerations  and  Momentum  Balances 

The  observations  discussed  in  the  previous  sections  depict  vigor- 
ously time-dependent  fluctuations  in  the  equatorial  current  system. 

To  better  understand  the  dynamical  processes  involved,  estimates  of 
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terms  in  the  equations  of  motion  have  been  made.  In  their  usual  form, 

2 

these  equations  are  in  units  of  acceleration  (cm/sec  ). 

Terms  which  have  been  examined  include: 


u, , vu,t  B yv,  v , 6 yu,  v u , etc. 

t’  y 3 ’ t ’ 3 V 7.7. 


The  accelerations,  like  the  velocities  from  which  they  are  derived, 
are  quasi-per iodic , so  that  the  instantaneous  values  are  considerably 
greater  than  mean  values. 

Typical  magnitudes  for  the  local  and  advective  time  derivatives 
-4  2 

of  u and  v are  10  cm/sec  These  are  in  agreement  with  the  values 
found  by  Knauss  (1966)  for  the  advective  derivatives.  In  the  region 

just  above  the  core  of  the  Undercurrent,  dissipative  terms  of  the  same 

...  2 
magnitude  are  encountered,  assuming  vertical  eddy  viscosity  of  50  cm/ sec 

The  choice  of  a value  of  vertical  viscosity  is  difficult:  Jones  (1973) 

reports  a strong  dependence  of  the  viscosity  on  depth.  In  the  light 

of  our  observations,  viscosity  is  likely  to  be  time-dependent  as  well. 

Poleward  of  0.75°  latitude,  the  Coriolis  acceleration  becomes  important. 

It  is  not  practicable  to  estimate  zonally  dependent  accelerations 

such  as  uu^  with  the  available  data,  but  they  are  likely  to  be  important. 

- 5 2 

Katz  (1976)  reports  a mean  value  of  p of  about  10  cm/sec  in  acceler- 

x 

ation  units.  Large  variances  of  p indicate  that  instantaneous  values 

x 

are  considerably  higher.  [Weisberg  (personal  communication)  reports 
meridional  pressure  gradients,  based  on  STD  data,  which  vary  on  similar 
time  scales.  ] 

2 

A windstress  of  the  order  of  0.5  dyn/cm  would  produce  the  observed 
value  of  v^  of  10  ^ cm/sec^  (assuming  a t /H  stress  model).  It  is 
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quite  probable  that  the  magnitude  of  the  time-dependent  windstress 

exceeded  this  value  at  R/V  ISELIN  during  Phase  II  of  GATE  (Fig.  4). 

Vertical  advection  terms  appear  to  be  important  in  regions  of  high 

vertical  shear  of  horizontal  velocity.  In  particular,  at  a depth  of 

-4 

50  in,  where  u reaches  a maximum,  wu  reaches  a magnitude  of  10 
£ z 

cm/ sec^ . 

Many  terms  in  the  momentum  equations  appear  to  be  of  comparable 
magnitude  near  the  equator,  and  it  is  not  clear  whether  there  are  any 
principal  balances.  Consequently,  an  attempt  has  been  made  to  systema- 
tically test  certain  balances.  Clearly,  since  we  cannot  evaluate  all 
terms  in  the  equations,  ail  possible  balances  cannot  be  considered. 

If,  however,  several  terms  seem  to  cancel  each  other  consistently,  the 
importance  of  the  terms  in  a dynamical  process  will  have  been  established. 
As  it  turns  out,  one  firm  conclusion  can  be  drawn  from  the  calculations 
to  be  described. 

We  test  a balance  among  terms:  A.(r,  t),  i = 1,  N 

where:  r = position  vector 

t = time 

A.  = term  in  equation 

l M 

by  computing  a coefficient: 

r ? "1^  ~ 2 k 

o(r ) = 1 - (I  A ) j / E (A.  ) * 

G l 1 J 11 

where  .j, 

(_)  ; | ( ) dt 

o 

and  o is  confined  to  the  interval  (0,  1).  A value  of  1 implies  a 
perfect  balance,  while  a value  of  0 results  for  maximum  imbalance. 
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We  are  calculating  the  complement  of  the  normalized  root  mean  square 
(RMS)  residual  of  a sum  of  quantities  which  are  expected  to  cancel  to 
a degree.  Balance  coefficients  were  computed  for  sets  of  terms  such 
a s : 


a) 

u 

+ 

Vll 

3vv 

t 

y 

b) 

u 

+ 

V'J 

Byv  - v u 

t 

y 

V 

c) 

V 

t 

+ 

3yu 

etc . 

Clearly,  any  calculations  involving  v ^ must  be  treated  with  caution, 
as  its  value  is  questionable. 

The  primary  result  shows  that,  at  the  depth  of  the  core  of  the 

/ 

Undercurrent,  particularly  at  latitudes  of  maximum  meridional  shear, 
there  is  a dominant  contribution  to  the  local  acceleration  by  the  meri- 
dional advection  of  zonal  momentum.  This  balance  coefficient  (computed 
for  a), above)  is  depicted  in  Fig.  18.  Equatorward  of  0.75°  latitude, 
the  Coriolis  term  is  of  secondary  importance . It  is  emphasized  that 
this  calculation  is  more  than  just  an  estimate  of  the  magnitude  of 
the  nonlinear  term,  vu  . It  represents  the  systematic  contr ibut ion 

y 

to  the  local  acceleration  by  that  term.  Other  balances,  including  those 
for  the  meridional  equations,  indicate  likewise  that  non-linear  terms 
are  important  in  most  of  the  latitude-depth  region  spanned  by  the  data. 
Unfortunately,  while  inclusion  of  more  terms  tends  to  improve  these 
balances,  they  by  and  large  do  not  rise  above  60%,  thus  limiting  the 
interpretation  of  the  results. 

The  conclusion  to  be  drawn  from  these  calculations  is  that  mean 
flow  must  be  considered  when  attempting  to  model  the  observed  processes. 
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At  the  very  least,  meridional  advection  of  zonal  momentum  and  of  the 
associated  geostrophic  meridional  pressure  gradient  must  be  included. 

2.3  Vertical  Motions  and  Meridional  Circulation 

An  important  aspect  of  equatorial  current  observations  is  the  meri- 
dional circulation,  i.e.,  north-south  and  vertical  flow.  Most  discussions 
of  these  flows,  however,  are  associated  with  steady-state  models  (e.g., 
McKee,  1973).  In  such  models,  a secondary  circulation  is  required  t:o 
balance  the  Ekman  divergence  resulting  from  the  zonal  wind  forcing. 

The  balance  is  effected  through  dissipative  terms  which  depend  on  the 

shear  in  the  meridional  flow.  Vertical  velocities  in  McKee's  model 

-2 

are  of  the  order  of  10  cm/ sec.  The  model  is  for  a homogeneous  ocean, 
however,  and  in  models  which  include  stratification,  the  magnitudes 
of  the  vertical  motions  and  meridional  flow  are  considerably  lower. 

Figures  6 to  10  represent  instantaneous  values  of  the  meridional 
circulation  derived  from  the  digitized  y-t  maps.  The  w-component  was 
computed  as  follows: 

z 

r> 

w(y,z,t)  = - J VI(.  vH  dz 
o 


where:  w( z = 0)  » 0 
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If  is  considered  important,  however,  we  must  make  some  assumption 
about  the  longitudinal  structure  based  on  the  observations.  Our  best 
g..<'ss  is  for  a we st ward- pr opagat ing  disturbance  with  a phase  speed 

of  1 hO  cm/ sec . 

If  we  let 

u=uk(x+ct)]  =u(  | ), 

/ 

the*  U s KCU 

t 

/ 

u * ku 
x 

/ _ 

u u - 


Ir.i  lusion  of  thr  x-btructure  results  in  values  of  w which  are 
smaller  bv  about  SO*  than  those  computed  from  v alone.  However,  the 

y 

overall  v-z  structure  remains  qua  1 i ta t ive ly  the  same.  Typical  magni- 

_o 

tudes  of  w were  found  to  be  10  A cra/sec  and  greater.  Errors  are  such 
that  w can  only  be  estimated  to  within  a factor  of  about  two.  The  main 
purpose  of  these  calculations  is  to  show  the  y-z  structure,  or  the 
meridional  flow,  and  to  illustrate  the  intense  time  dependence. 

There  is  a systematic  structure  to  the  derived  meridional  circu- 
lation in  the  y-z  plane,  and  it  is  quasi-per iodic  in  time.  Due  to  large 
time-dependent  amplitudes  and  a short  record  length,  the  time-averaged 
value  of  this  circulation  is  not  significantly  different  from  zero. 

Thus,  comparisons  with  meridional  circulations  implied  in  steady-state 
models  are  not  meaningful.  On  the  other  hand,  the  oscillatory  nature 
of  w and  v (and  well  as  u)  are  consistent  with  a wavelike  phenomenon. 

The  meridional  circulation  is  essentially  consistent  with  the  variation 
in  the  position  and  with  the  shape  of  the  zonal  current  field.  The 
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results  imply  that  the  time-varying  meridional  circulation  is  part  of 
the  same  wavelike  phenomenon  that  is  manifest  in  the  lateral  oscillation 
of  the  Undercurrent.  Based  on  observational  results  alone,  it  is  not 
clear  whether  the  oscillations  are  solely  meanders  in  the  Undercurrent 
or  are  manifestations  of  some  wave  phenomenon  which  extends  beyond  the 
north-south  limits  of  the  Undercurrent 

Figure  19  is  a plot  of  time  series  of  meridional  transports  in 
the  two  previously  defined  layers,  and  their  sum: 


M1  = H1V! 
M2  = »2V2 


M„  = M,  + M„ 
T 1 2 


The  total  (M^)  shows  the  16-day  periodicity  even  more  clearly 
than  in  the  Undercurrent  or  surface  flow.  corresponds  roughly  to 

the  external  mode  of  a linear,  two-layer  system. 


2 .4  Wind  Observations 

In  Figs.  2 to  4,  similarities  in  scales  of  variability  can  be  seen 
between  wind  and  current.  Correlation  coefficients  for  the  v-components , 
calculated  for  several  time  lags,  are  presented  in  Fig.  20a.  A lag  of 
one  day  yields  a high  degree  of  correlation  between  about  1°S  and  0.75°N 
latitude.  Correlations  between  the  zonal  components  are  much  lower, 
except  near  1°N,  for  a time  lag  of  5 to  6 days.  Although  these  results 
cannot  be  considered  to  be  statistically  significant,  they  suggest  that 
a high  correlation  between  meridional  components  of  wind  and  current 
is  plausible. 
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Time  series  of  wind  observations  from  two  Soviet  ships  at  the 
equator  have  been  examined.  The  ships  were  R/V  AK.  KURCHATOV  and  R/V 
PASSAT.  Observations  were  obtained  at  23.5°W  and  10°W,  respectively. 
Figure  20  b is  a typical  periodogram  computed  from  the  PASSAT  data. 

A significant  amount  of  energy  exists  in  the  v-component  at  periods 
ranging  from  5 days  to  several  weeks.  Gaps  in  the  records  between 
phases  necessitated  separate  analysis  of  each  phase,  thus  reducing  the 
quality  of  the  results. 

Krishnamurti  (1975b)  inferred  the  wind  field  at  850  mb  from  satel- 
lite photographs  of  clouds  during  the  entire  GATE  period.  Subsequent 
analysis  of  the  winds  revealed  simultaneous  surges  in  the  Tradewinds, 
with  a dominant  period  of  14  to  15  days  and  longitudinal  wavelength 
of  the  order  of  2500  km.  At  the  equator,  these  surges  were  primarily 
from  the  southeast,  with  a maximum  near  25°W.  Amplitudes  of  the  surges 
were  about  3 m/ sec , with  maximum  intensities  up  to  8 m/sec.  These 
values  are  in  agreement  with  those  observed  from  R/V  ISELIN  (recall 
Fig.  4).  Krishnamurti  (1975a)  cites  other  authors  reporting  2 to  3- 
week  oscillations  in  tropical  wind  systems  (Wallace  and  Chang,  1969; 
Yanai  and  Murkami,  1970;  Kousky  and  Wallace,  1971;  Krishnamurti  et  a 1 , 
1972;  Gruber,  1974). 

In  summary,  wind  observations  during  GATE  imply  the  existence  of 
time  variability  of  the  meridional  wind  component  at  periods  of  2 to 
3 weeks,  with  longitudinal  scales  corresponding  to  wavelengths  of  2000 
to  3000  km,  A very  simple  argument  will  now  be  presented  to  show  that 
a modest  windstress  applied  on  this  time  scale  can  produce  an  accelera- 
tion resulting  in  a velocity  amplitude  comparable  to  that  observed 
during  the  PCM  experiment. 
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Assume:  t = 0.25  dvn/cm.  Let  the.  windstress  act  as  a body  force 

over  a layer  of  depth,  H = 40  m.  A period  of  16  days  has  a quarter- 
cycle  of  4 days  (the  time  during  which  the  response  should  reach  maximum 
amplitude).  Thus, 

4 days  , , 

<*  (y) 

v = — — at  = 22  cm/' sec.  (o  55  1) 

t = 0 0H 

This  value  is  indeed  of  the  same  order  as  the  observed  amplitude. 

The  stress  model  described  is  th“  one  used  in  Chapter  3 in  the  dynamical 
model  to  be  considered. 

2 . 5 Supplementary  Observations 

At  23.5°W  and  10°W,  north-south  arrays  of  Soviet  BPV-2  (or  Alekseev) 
moorings  were  deployed  by  R/V  KURCHATOV,  R/V  PASSAT,  R/V  VIZE,  and  R/V 
ZUBOV.  These  moorings  were  maintained  during  all  three  phases  of  GATE. 
Each  mooring  carried  about  10  BPV-2  current  meters  at  depths  ranging 
from  15  to  700  meters.  They  sampled  current  speed  and  direction  every 
half-hour.  In  a preliminary  GATE  report  and  in  TROPEKS  *74,  Bubnov 
(1975  a & b)  analyzed  the  observations  from  23.5°W.  The  resulting 
la t i tude-t ime  maps  of  velocity  components  show  the  same  oscillatory 
behavior  in  the  equatorial  current  system  as  observed  from  R/V  ISELIN 
and  R/V  A.  V.  HUMBOLDT  (Fig.  21).  In  a later  study,  Bubnov  (1976) 
presents  the  results  of  a spectral  analysis  done  on  selected  time  series 
from  both  23.5°W  and  10°W.  The  principal  result  was  the  existence  of 
a pronounced  energy  peak  for  a period  of  16  to  18  days  at  and  near  the 
equator  at  depths  up  to  100  m.  There  was  also  evidence  of  higher  fre- 
quencies, including  semi-diurnal  tides.  Off  the  equator,  at  10°N  and 
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23.5°W,  a dominant  peak  was  found  at  the  local  inertial  period  of  68 
hours.  Bubnov  concludes  that  the  several-week  period  corresponds  to 
an  equatorial  inertial  oscillation.  It  is  noted  that  the  inertial 
period  at  one  baroclinic  radius  of  deformation  (100  to  300  km)  from 
the  equator  is  of  the  order  of  2 to  3 weeks. 

Meincke  (1975)  obtained  results  from  three  current  meter  moorings 
at  29  W at  latitudes  of  0°1'S,  0°49.5'S,  and  0°25'N.  The  mooring  at 
0°25'N  yielded  one  record  at  a depth  of  18  m.  At  0°1'S,  records  at 
six  levels  were  obtained:  18,  75,  96,  127,  308,  and  499  meters.  At 

0°49.5'S,  records  at  five  levels  were  obtained:  75,  96,  127,  308,  and 

499  meters.  The  instruments  operated  throughout  Phase  II.  Spectral 
analysis  yielded  dominant  energy  peaks  at  semi-diurnal  periods.  Less 
pronounced  peaks  at  24.8  hours  were  also  present.  Records  were  not 
of  sufficient  duration  to  yield  spectral  estimates  for  periods  longer 
than  about  8 days,  but  examination  of  plots  of  the  v-component  reveals 
a pronounced,  quasi-per iodic , 16-day  signal  with  an  amplitude  of  about 
40  cm/sec  to  a depth  in  excess  of  100  m.  Tidal  signal  amplitudes  were 
typically  10  to  15  cm/ sec . 

Krishnamurti  (1976)  reported  on  the  objective  analysis  of  sea 
surface  temperatures  (SST's)  based  on  satellite  data  and  ship  obser- 
vations. Spectral  estimates  based  on  the  data  from  the  entire  GATE 
period  show  fluctuations  with  periods  of  3,  7,  16,  20,  and  38  days, 
as  well  as  higher  frequencies,  including  tidal.  East-west  scales  of 
several  thousand  kilometers  are  evident  in  the  analyzed  maps  of  SST. 

Brown  (1977)  reports  similar  spatial  scales  in  satellite-derived 
SST  data  from  the  equatorial  area  during  GATE. 
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Current  profiles  tanen  by  Hallock  (1975)  aboard  R/V  KURCHATOV  at 
the  equator  and  at  23.5°W  show  evidence  of  tidal  frequencies  (the  samp- 
ling interval  was  6 hours)  and  several-week  lateral  and  vertical  fluc- 
tuations in  the  position  of  tl  - Undercurrent.  The  duration  of  this 


series  of  profiles  was  about  two  weeks. 

Ail  these  observations  imply  that  the  16-day  oscillation  is  a real 
and  prominent  feature  of  the  equatorial  current  system,  and,  while 
signals  of  higher  frequency  are  oresent,  the  PCM  data  do  not  seem  to 
be  severely  alaised  bv  these  signals. 


2 . 6 Summary  and  interpretation 

The  GATE  PCM  observations  offer  strong  evidence  of  north-south 
oscillatory  motions  with  a 2 to  3-week  period  in  the  equatorial  Atlantic. 
Associated  with  these  fluctuations  are  vertical  motions,  as  inferred 
from  calculations  of  vertical  velocitv  which  are  also  periodic  in  time. 

The  meridional  circulation  is  consistent  with  changes  in  the  latitude- 
depth  structure  of  the  zonal  current.  The  existence  of  a 16-day  oscil- 
lation is  corroborated  by  independent  measurements  at  longitudes  ranging 
o o 

from  10  W to  29  W during  all  three  phases  of  GATE.  The  principal  features 
of  these  fluctuations  are  retained  when  the  data  are  averaged  into  two 
layers,  the  interface  between  the  layers  corresponding  to  the  isotach 
of  zero  zonal  velocitv. 

Dynamical  calculations  revealed  that  advective  and  frictional  terms 
in  the  equations  of  motion  are.  probably  important  near  the  equator, 
and  in  particular,  meridional  advection  of  zonal  momentum  is  pronounced 
in  the  level  of  the  Undercurrent. 
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The  role  of  frictional  processes  is  found  to  be  important,  par- 
ticularly in  regions  of  large  vertical  shear  of  horizontal  flow.  It 
is  difficult,  however,  to  quantitatively  determine  its  effect  in  the 
present  case  due  to  the  uncertainty  of  a value  for  vertical  eddy  vis- 
cosity. 

It  is  shown  that  the  assumption  of  a reasonable  windstress  acting 
over  the  upper  of  the  two  defined  layers  is  capable  of  producing  the 
observed  northward  speeds.  Furthermore,  both  the  frequency  content 
of  the  wind  field  and  its  correlation  with  the  surface  flow  suggest 
wind  forcing  on  the  scales  of  the  observed  16-day  fluctuation. 

Excitation  by  a time-varying  wind,  with  possible  strong  interac- 
tions with  mean  zonal  currents,  constitutes  a plausible  mechanism  for 
generation  of  the  several-week  oscillation.  This  conclusion  provides 
the  motivation  for  the  theoretical  discussions  and  numerical  experiments 


to  follow. 
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3 . Theoretical  Considerations and  Nume ricai  Experiments 


3 . 1 Introduction 

As  discussed  in  Chapter  2,  several  investigations  of  the  response  of 
the  equatorial  ocean  to  winds  have  been  carried  out  using  a two-layer  ap- 
proximation for  the  vertical  structure.  The  primary  reason  for  choosing 
such  a system  is  that  it  is  the  simplest  one  which  can  accomodate  the 
gross  features  of  the  structure  of  the  equatorial  Currents.  A single 
layer  would  necessitate  extending  the  Undercurrent  to  the  surface  every- 
where and  would  not  allow  for  the  simulation  of  some  of  the  features  in 
the  observations.  For  example,  the  southward  shift  of  the  core  of  the 
Undercurrent  four  days  after  the  onset  of  the  northward  flow  at  the  sur- 
face would  not  be  possible  in  a single  layer  system.  The  two-layer 
system  also  allows  severs'  options  for  setting  the  density  structure 
and  for  the  application  f forcing.  As  will  be  shown  later,  the  layered 
system  in  many  ways  resembles  a continuously  stratified  system,  such  as 
the  one  treated  by  Moore  (1968). 

The  system  to  be  considered  may  be  referred  to  a " 2 layer  model. 

It  consists  of  two  active  layers  of  finite  mean  thickness  overlying  an 
infinitely  deep  or  quiescent  layer  (Fig.  22).  Each  layer  is  treated  as 
a shallow  fluid  of  constant  density.  Densities  increase  with  mean  layer 
depth.  A property  of  the  bottom  layer  thus  defined  is  that  horizontal 
pressure  gradients  and  horizontal  motions  are  everywhere  vanishingly 


small.  The  elevation  of  the  lower  interface,  <n  , is,  therefore,  com- 

3 

pletely  determined  by  the  elevations  of  the  other  two.  The  horizontal 
extent  of  the  system  is  from  15°N  to  15°S.  For  the  class  of  motions 


considered,  and  duo  to  the  fact  that  the  2'^- layer  configuration  has  no 
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barotropic  mode,  boundaries  at  these  latitudes  have  virtually  no  effect 
on  the  short-term  response  near  the  equator.  The  east-west  structure  is 
approximated  by  simple  harmonic  dependence. 

A mean,  geostrophically  balanced  zonal  current  is  included  in  the 
model.  It  consists  of  a westward  flow  in  the  surface  layer  and  an  east- 
ward flow  in  the  Undercurrent  layer.  Meridional  dependence  of  the 

2 

currents  is  Gaussian  (i.e.,  e-y  ) with  an  e-folding  distance  of  300  km 
in  the  surface  layer  and  100  km  below.  The  current  fields  thus  defined 
are  fixed  in  time  and  independent  of  longitude.  They  enter  the  equations 
as  latitude-variable  coefficients  simulating  some  of  the  non-linear 
terms.  Due  to  these  properties  it  is  possible  to  Fourier  transform  the 
equations  in  longitude  and  consider  the  problem  for  a single  longitudinal 
wave  number  at  a time. 

Forcing  is  applied  by  assuming  an  acceleration  in  the  form  of  a body 
force,  t/H,  where  T is  the  windstress  and  H is  the  layer  over  which  it 
acts.  The  windstress  in  this  study  is  a simple  function  of  time  and  is 
independent  of  latitude,  except  near  the  boundaries.  Longitudinal  de- 
pendence is  implicit  in  the  wave  number  of  the  harmonic  component  con- 
sidered. Attention  has  been  restricted  to  a meridionally  forced  problem 
since  it  is  the  simplest  form  of  forcing  that  will  produce  a latitudi- 
nally  symmetric  v-component,  such  as  that  observed,  in  the  response.  To 
obtain  such  a response  with  only  a zonal  forcing  component  would  require 
a nonzero  windstress  curl.  Clearly,  the  observed  wind  field  is  much 
more  complex  than  the  assumed  forcing  function,  but  the  aim  here  is  to 
seek  a possible  mechanism  for  wave  generation,  not  to  explore  all 
possible  cases. 

In  the  following  sections,  the  equations  for  this  system  are 
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presented  and  put  Into  a form  to  be  integrated  numerically.  Also,  the 
linear  case  has  been  examined  analytically  for  comparison  with  the  obser- 
vations and  to  check  numerical  results.  Finally,  the  numerical  simula- 
tion itself  is  discussed  and  cases  with  mean  flow  are  compared  with  other 
results . 

3 . 2 Derivation  of  Model  Equations 

I begin  with  the  linearized,  inviscid,  shallow-fluid  equations  on 
an  equatorial  beta-plane  (i.e.,  the  Coriolis  parameter  is  directly  pro- 
portional to  the  latitude)  for  a two-layer,  divergence-coupled  system 
overlying  a third,  infinitely  deep,  motionless  layer.  This  system  is 
essentially  the  same  as  that  used  by  Philander  (1976).  Included  is  a 
f ixr  ' , geostrophically  balanced  mean  zonal  current  in  each  of  the  active 
layers.  I have  assumed  that  the  non-linear  terms  are  small,  except 
where  interaction  with  the  idealized,  longitudinally  independent  zonal 
currents  are  concerned.  Thus,  for  the  fluctuations  (dependent  variables), 
the  Rossby  number  is  small.  A Rossby  number  is  defined,  however,  which 
relates  to  the  magnitude  of  the  mean  flow.  Wind  forcing  is  introduced 
as  a body  force  which  acts  uniformly  over  the  upper  layer. 

For  convenience,  the  following  constants  are  defined: 


On  the  equatorial  beta-plane,  the  Coriolis  parameter  takes  the  form: 

f = 20  sin  5 = —f  s By 

R 
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where 


also : 


n = 
§ = 
R = 

y = 

X = 


(W  = 


8 = 


(F,  G)  = 


rotation  rate  of  earth 
latitude  (angle) 
radius  of  earth 
latitude  (distance) 
longitude  (distance) 
horizontal  velocity  in  layer  i 
acceleration  of  gravity 
upper  layer  forcing  function 


Following  Fig.  22  for  other  variable  definitions,  the  equations 
of  motion  for  this  system  are: 


U1  + Vi  + Vi 

- ByVl  = 

t X 

y 

V1  + ulvl  + 0yu1 

= -gT^  + G 

1 t X 

y 

@yui  = -g^ 

y 

+ - \ )+  v 

<\  - \ >+ 

X X 

y y 

u2  + G2U2  + v2G2 

- Byv  = -g 

t x y 

°2 

V2  + U2V2  + Syu2 

I^(Oi  ^ 

t X 

o2  y 

Byu 

= =&  (pi\ 

2 

°2 

+ F 


L + Vi  ) = ° 

x y 


Eqs.  (1) 
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— T TL  + a\  + u~  (~ V ^1  + <*11.  ) + v„  f— ir  11  + oft  „ ) + 

pi  2 2 p 1 2 2 D 1 2 

t t x x y y 

(h9  + ^ t 2 ) (u?  + v2  ) = 0 


The  boundary  conditions  are  that  the  meridional  component  of 
velocity  vanish  at  the  zonal  boundaries  at  15°  latitude,  and  all 
variables  are  taken  to  be  periodic  in  the  selected  east-west  wave 
number.  The  system  is  initially  at  rest  (except  for  the  mean  flows) 
and  is  integrated  in  time  with  an  applied  forcing  function.  Hie 
boundary  conditions  are  expressed  as  follows: 


At  y = ± 15°  latitude 


u 


-g 


1 

x 


V 


1 


= 0 


u 


2 


t 


) 


u ^ (x  + \)  = u1  (x) , etc. 


ikx 

For  « e ; 


k 


iz 

\ 
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(The  mean  flow  is  prescribed  in  such  a way  as  to  be  virtually  zero 
at  and  near  the  zonal  boundaries.  Hence,  they  need  not  be  included 
here . ) 

Non-dimensional  variables  and  associated  scale  constants  are 
defined  as  follows: 

(u^ , u2,  V , v2)  su  (u^,  u9',  v',  v2') 

(X,  y,  \)  = L (x',  y',  \'),  H1/H2  5 H 
t -Tt',  ^ ^ SA2  ^ 

(u^  u2)  = U (u',  u2'  >»  \ 3 \ Pi  . * 2 S h \ 


Scale  factors  are  chosen  which  are  appropriate  for  the  class  of 
motions  being  investigated.  This  choice  is  largely  a matter  of 
convenience  since  no  terms  will  be  neglected  as  a result. 

Anticipating  accelerations  of  similar  magnitude,  but  not  neces- 
sarily of  the  same  phase  in  the  two  layers,  pressure  gradients  in  the 
two  layers  are  assumed  to  be  of  comparable  magnitude  as  well.  Hence, 
the  same  scale  factor  is  chosen,  appropriately  weighted  for  the 
density  differences,  for  the  two  interfaces,  Ti  and  ^ . 


b9  (A.  , A. ) 

(A2,  A2)  = (A , A ) h - -jt 
/ 

(A. , A ) = — — (A,  A) 

I I 0O 


or 
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An  equatorial  Rossby  number  is  defined  which  is  associated  with 
the  magnitude  of  the  mean  current: 


R 

o 


= ii/BL2 


In  the  momentum  equations,  relations  among  U,  A,  and  L are 
chosen  to  be  consistent  with  quasi-geostrophic  scaling: 


(U,  U) 


(A,  A ) 
BL" 


A baroclinic,  equatorial  radius  of  deformation  is  defined: 


R 


d 


Cg  = ^ g P/H1/p2 


Thus,  (A,  A)  = (U,  U)  H^BR2 


Finally,  defining: 


* — / / // 
ft  = Pj  0 / 020 


J 
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and  assuming 


P1/o2  « 1. 

and  choosing  an  "inertial"  time  scale: 
T = 1/BL 


the  scaled  momentum  equations  are  divided  by  0l  (U,  U)  and  the  continuity 
equations  by  A/T,  which  yields: 


(U1U1 


Vi 


) - yv. 


+ R u v 
oil 

X 


+ yUj 


y“i 


G 

ubl 


V\  - "^2  + Ro  L“l  (y\  - \ ) + v!  (y\  - \ )]  + 

t t xx  v v 


Eqs. (2) 


Rd 


+ Rq  (y\  - j (uT  + v ) - 0 

L x y 


u2  + Ro  (G2u2  + v2G2  ) - yv2  = 

t x y 


X X 


v2  + RoU2v2  + yu2  = -^l  - \ 

t x y y 


y“2  ■ *\  - \ 
y y 
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+ 0'712  + R0  ^u2  + a\  ) + v2  (6^  + ) 1 + 

t t xx  y y 

r 2 

H + Rq  (6^  + c*^)]  (u2  + v2  ) = 0 

L x y 

It  will  be  shown  later  that  taking  L = R.  is  a convenient  choice 

a 

since  R^  turns  out  to  be  the  equatorial  trapping  scale  for  the  motions 

involved.  L has  been  left  in  Equations  (2)  to  illustrate  how  the 

divergence  of  horizontal  velocity  enters  the  dynamics.  For  example, 

2 

for  T.  » Rj”,  the  motion  is  essentially  non-divergent . The  remaining 
free  parameters,  U,  U correspond  to  typical  amplitudes  of  the  fluctua- 
tion velocity  field  and  the  mean  flow,  respectively.  U is  set  equal 
to  for  convenience  since  it  appears  only  in  the  forcing  terms. 


Ia_G 

U3L 


F,  G F,  G 

qi  = 2 q 

(BCV  6 Rj 

g d 


All  other  constants  are  determined  by  the  static  structure  of  the 
system  (i.e.,  the  density  field)  and  the  rotation  rate  of  the  earth. 

From  the  GATE  data  as  well  as  previous  observations  (Stalcup  and 
Metcalf,  1966;  and  Rinkel,  1965),  layers  can  be  defined  such  that  the 
following  parameters  can  be  evaluated. 

H1  = 40m,  H2  = 80m,  - H = 2 

p'  = p"  = 0.002  - a = 2,  s =1 
P1  as  P2  fa  p3  » 1 


(aside  from  differences) 
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This  yields : 


Cg  = 88.5  cm/sec  = U 

L = R,  = 197  km 
d 

T =2.57  days 

The  time  scale  can  be  defined  as  the  period  divided  by  2 jt  . Thus, 

Per  = 2 ttT  = 16  days 

The  baroclinic  adjustment  time  scale,  T,  is  compatible  with  the 
dominant  time  scale  present  in  the  observations.  If  the  system  were 
treated  as  one  active  layer  120  m thick  with  an  overall  density 
difference  of  0.004,  one  finds: 

C = 217  cm/sec,  L = = 308  km,  T = 1.64  days.  Per  = 10.3  days 

8 

Thus,  a rather  large  change  in  problem  parameters  results  in  a relative- 
ly small  change  in  the  time  scale.  Since  the  choice  of  and  hence, 

C ,is  somewhat  arbitrary,  these  results  may  seem  fortuitous.  It  is 
8 

shown  in  Section  3.3;  however,  this  is  not  the  case. 

Rq  is  quite  dependent  on  problem  parameters  for  U in  the  range  of 
observed  values.  Table  2 illustrates  this  point. 


Table  2 : Rossby  No.  as  Function  of  Gravity  Wave  Speed  and  Current 
Velocity 


U (cm/sec) 

C (cm/sec) 

R 

o 

10 

88.5 

0.113 

10 

217.0 

0.046 

25 

88.5 

0.281 

25 

217.0 

0.115 
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Harmonic  dependence  in  longitude  (x)  is  now  introduced. 


u^(x>  y,  t)  - Re  ux  (y,  t)  e ,y  = u^  cos  kx  - u sin  kx 

r 


where 


k = 2tt/X,  X = east-west  wavelength 


Applying  scaling: 


X = LX'  , k = k'/  L 


Since  the  equations  are  mathematically  linear,  even  when  mean  flow  is 
included,  once  the  harmonic  dependence  in  x has  been  introduced,  they 
can  be  Fourier  transformed  in  that  coordinate  so  that  only  one  wave 
number  component  need  be  considered  at  once;  i.e.  : 


u + uu  = F 
t x 

00 

r-  ik.x  ik.x 

u e + ik  uue 

j=l  1 


where  k.  are  such  as  to  make  the  set 

ik  x 

X1 , where  k = mk  . Multiply  by  e m , integrate  over  X,,  and  divide 

ik  x m 1 1 

, m 

by  e : 

-*  u\  + ik  uu  = F 
t m 


u = ue 
1 = 1 


ik  .x 
1 


ik  .x 


Fe 


j=l 


ik  .x 


e > j = 1 > ® orthogonal  over 
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The  resulting  transformed  equations  are  (omitting  primes  and  ~ ) : 


- - 2 3 

u.  + R (iku.u  + v u ) - yv,  = -ikT  + F/8  R, 

l.ollll  1 1 d 

c y 

v + R iku.v  + yu  = -"H  + G/82R2 

1 o 1 1 1 1 d 

t y 


y\  - \ + RQ  iku^v^  - + V;L  (Y\  - \ ) ! + 

t t y y 


1 + Ro(Y\  - (iku^  + v1  ) = o 

y 


Eqs.  (3) 

U2  + Ro^ik^2U2  + V2^2  ^ ' yv2  = -ik(\  + "V 

t y 


v2  + Ro  ikQ2V2  + yu2  = \ 

t y y 


5^  + «'n2  + Rq  iku1  (^'n  + a'Tl2)  + v2  (gTL  + Of^  ) + 

t t y y 


H + R (git  + a\)  (iku  + v ) = 0 
o 1 Z Z Z 

y 


(The  mean  flow  equations  are  unchanged.  In  fact,  they  simply  represent 
the  relationships  of  U to  T|  and  do  not  enter  explicitly  into  the  solution 
of  the  problem;  i.e.,  the  mean  flow  is  prescribed.) 


These  six  equations,  with  the  previously  stated  boundary  conditions, 
constitute  the  system  to  be  integrated  numerically  in  time.  Prior  to 
this,  the  integration,  an  analysis  of  the  linear  equations  for  free  waves 
(i.e.,  for  Ro=0  and  no  forcing)  will  be  carried  out  to  ellucidate  numeri- 


cal results. 
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3.3  Solution  of  the  Linear  Problem 

The  problem  of  linear,  equatorial  waves  has  been  dealt  with  at 
length  by  Lighthill  (1969),  Matsuno  (1966),  Moore  (1968),  and  others. 
The  purpose  of  Section  C is  to  put  the  previously  derived  system  of 
equations  into  a form  which  is  compatible  with  existing  theory  for  the 
case  where  Ro«l.  The  results  of  the  linear  analysis  are  compared 
with  the  observations  and  also  with  the  numerical  results  for  a zero 
mean  flow  as  a check  on  the  validity  of  the  model. 


3.3.1  Derivation  of  the  ver tical  normal  modes 

The  objective  of  the  following  procedure  is  to  take  the  six 
coupled  Equations  (2),  with  Ro=0  for  the  two-layer  system,  and  define 
two  sets  of  normal  (uncoupled)  equations  which  can  be  solved  separately.* 
Each  of  these  sets  is  equivalent  in  form  to  the  other  and  to  the  equa- 
tions of  a barotropic  (one-layer)  system.  The  subsequent  analysis  is 
identical  to  that  used  by  Matsuno  (1966)  yielding  a dispersion  rela- 
tion for  each  mode. 

I begin  with  the  redimensionalized  equations  where  the  mean  flow 
terms  have  been  omitted: 


u - SyvL  = -g^  + F 

*t  x 


V + Syu  = -gT^  + G 
t 1 y 


\ - \ + ht  (U]  + Vl  ) = 0 

t t x y 


By.  . -g(-l  T, 

2 x 2 x 


* The  procedure  used  is  that  outlined  by  Lighthill  (1969) 


Boundary  conditions  are  taken  at  y = ± oo  and  require  all  fields  to 

vanish.  Following  Lighthill  (1969),  thickness  variables  are  defined: 

hl  = \ ‘ V h2  = \ " "S’  h3  = ^ 


(Constants  such  as  H have  been  omitted  from  h since  the  thicknesses 
enter  only  as  derivatives.) 

As  previously: 


h3  " 


' Dlhl  '°2h2 


resulting  with: 


u,  - Byv  = ^ ' o (1  - ) h + p (1  - ) h + F 

Ltt  Xt  °l  u 1 °3  Xxt  1 p3  2xtJ 


u?  ■ 9yv2  = 7s  pi  (1  - -r  > h!  + p,  a - ~ > h + g 

2tt  2t  °2  " 1 p3  Xxt  2 p3  2xt  c 


I 


43 


hl  = "H1  (U1  + V1  > 

xt  xx  xy 


h2  = -H2  (u2  + v2  ) 

xt  xx  yx 


where  the  momentum  equations  have  been  differentiated  with  respect 
to  time  and  the  continuity  equations  with  respect  to  x (the  procedure 
for  the  v-momentum  equations  is  completely  analogous). 

Defining  a matrix: 


p3 


p3  ‘ 


A = 


a . . 

n 


( 


3 


°2> 


' 0 o 


the  abo^e  equations  are  rewritten 


Jt. 


- 9y 


v . = - g / a . . h . + F . . = g)  a . . H . (u.  +v. 

jt  ji  i . jt  ji  i i i 

i=l  xt  J 1=1  J xx 


+ F. 

*y  Jt 


for  each  layer,  j. 


r 
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^ ~ ic 

Normal  mode  variables,  u , v , as  well  as  an  "effective"  or 
_ 

"equivalent"  deDth,  H , are  defined  such  that  the  following  is  true: 


-k  _ -k  -k  -k  -k  . -k 
u.  . - Byv  = -gH  (u  + v ) + F. 
tt  xx  xy  t 


Eq.  (4) 


The  layer  variables  are  defined  in  terms  of  the  normal  mode  variables 
by: 


( 


u,,  v ) = } Ck  (Gk,  vk)  , (F,  G)  = ) Ck  (f\  Gk) 

J J k=l  J k=l 


Since  the  normal  modes  are  defined  to  be  independent,  the  layer 
equations  can  be  written  separately  for  each  mode;  i.e.,  in  Eq . (4  ), 
substituting  for  u , for  mode  k: 


(ur  v5)  = cj  (“k’*k)  ^ (Fk>  ck)  = cj  (^k’  sk) 


This  results  with: 


2 

V 


,k.-k  Q -k.  V „ „k  ,-k 

. (u„.  - Byv  ) = -g  ) a . . H . C . (u 
i tt  t ./...  11  1 i xx 

J 1=1  J 


k '-k  + Vk  ) + Fk  Ck 

xy  t j 


-*bV  + ;k  > ck  + pk  ck 

xx  xy  j t j 


A new  matrix,  B,  is  defined: 


B = b . . = a . . II , 
j i J t i 


L 
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To  satisfy  the  above  equation,  the  following  is  required: 


> b..ck 
iSi  J1  1 


-k 

This  is  an  eigenvalue  problem  where  the  equivalent  depths,  H , are 

k 

the  eigenvalues  and  the  proportionality  constants,  , form  the  eigen- 
vectors of  the  matrix,  B.  Consequently,  for  each  vertical  mode,  k,  we 
haveH,  , j = 1,  2 and  a set  of  equations: 


utfc  ' = “8H  (uxx  + vxy) 


+ F 


vtt  + py~it  = -gH 


(u  + v ) + G 
xy  yy  t 


These  are  the  equations  of  a single-layer  or  barotropic, 
shallow  fluid,  such  as  that  treated  by  Matsuno  (1966),  whose  solution 
is  outlined  later  in  this  section.  The  above  procedure  is  applicable 
to  a system  with  any  number  of  layers,  as  described  by  Lighthill  (1969). 
The  case  of  continuous  stratification,  treated  by  Lighthill  (1969) 
and  Moore  (1968),  has  the  same  end  result,  namely,  the  above-normal 
mode  equations  and  an  equivalent  depth  which  is  proportional  to  the 

II  II  II 

Brunt-Vai sala  frequency.  The  matrix  problem  is  replaced  by  a Sturm- 
Liouville  problem  involving  the  solution  of  the  separated,  vertical 
structure  equation  for  the  system. 
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3-3.2  Calculation  of  normal  mode  parameters 

In  the  previous  subsection  we  derived  the  equation: 

BC  = HC 


where  C is  the  vector,  C j = 1,2  for  each  normal  mode,  k,  where 

J » 

the  superscript  has  been  omitted.  The  eigenvalues  are  found  by 
solving  the  following  equation: 


Det  [B  - Hi]  = 0 

where  I is  the  identity  matrix.  This  expression  results  in  the 
secular  equation: 

-2 

H - pH  + q = 0 


where 


P = b22  + bn 


q bub22  b21b12 


We  then  get: 


H =(P  + J 


P - 4q  )/  2 


Choosing  the  following  system  parameters; 


P3  ' d2  = °2  ‘ P1  = °-002 


Hl  = 40  m,  H2  = 80  m 
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and  assume:  p^  » ~ **  1 


(aside  from  differences) 


B = 


16 


16 


16 


H = (27.3,  4.69)  cm;  (mode  1,  mode  2) 


/- 

r _ 

1 

c = . 

.707 

- 

-.707 


Hereafter,  the  mode  corresponding  to  H = 27.3  cm  is  referred  to 
as  the  external  mode  and  the  other  as  the  internal  mode.  The  external 
mode  consists  of  a weighted  sum  and  the  internal  mode  to  a weighted 
difference  of- the  layer  velocity  components  in  the  light  of  the 
corresponding  eigenvectors,  C.  Both  modes  are  baroclinic  in  nature 
due  to  the  underlying  motionless  layer. 


u1  = (u^  + 2 u2)/  2 


external  mode 


«2  = (ux  * J 2 u2)/  2 


internal  mode 


The  external  mode  is  nearly  equivalent  to  the  single  mode  of  a 


J 
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system  consisting  of  one  active  layer  overlying  an  infinitely  deep, 
motionless  layer.  The  equivalent  depth  in  this  case  would  be  given 
approximately  by  (Lighthill,  1969): 

H - — H 


where  H is  the  active  layer  depth,  p is  the  density  difference,  and 
is  the  average  density  of  the  system. 

Associated  with  each  equivalent  depth  is  a gravity  wave  speed. 


C = gH.  This  is  not  to  be  confused  with  the  phase  speed  of  waves 

g J 

in  the  solution  (C^  = cr /k) , but  is  related  to  this  and  other  quanti- 
ties, as  will  be  shown. 

In  a later  subsection,  the  solution  to  the  horizontal,  normaj 
mode  equations  is  discussed.  Periods  and  phase  speeds  are  proportional 
to  1 /Jc  , Cg,  respectively.  It  is  useful,  therefore,  to  examine  the 
sensitivity  of  to  small  changes  in  problem  parameters.  This  has 
been  done  for  each  of  the  two  modes.  Values  of  all  pertinent  quanti- 
ties were  computed  for  a wide  range  of  layer  depths  and  density 
differences . 

From  these  calculations,  a subset  was  selected  (see  Table  3 
below)  as  being  representative  of  a reasonable  range  of  values. 

Table  3:  Range  of  Parameters  for  Error  Estimates 
H|=  (20,40,60)  m 
H2=  (40,80, 120)m 
— = (2,3,4)  X10-3 


= (1,2,3)  X10'3 
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From  these,  a range  of  values  for  H,  Cg,  and  2n/,/P Cg  was  derived 
for  each  mode.  The  results  are  summarized  in  Table  4 below. 

Table  4:  Range  ol  Normal  diode  Parameters 
External  Internal 


H = 14.3 

- 43.1 

OC 

CO 

- 9.53 

cm 

Cg  = 119 

- 205 

56 

- 97 

cm/ sec 

9 TT 

— — = 13.4 
/rc~ 

• g 

- 10.6 

20.3 

- 15.5 

days 

2-*  I 

It  turns  out  that  T = : 5 Cg  is  the  period  of  the  gravest 

horizontal  mode  (excluding  the  Kelvin  wave)  for  east-west  wave  number, 
k = 0.  This  calculation  shows  that  Cg  and,  to  a greater  extent,  T 
are  relatively  insensitive  to  realistic  variations  in  system  parameters. 

_i 

To  be  more  concise,  T varies  as  H , where  H depends  directly  on  layer 
depths  and  density  differences. 

According  to  Moore  and  Philander  (1975),  the  equivalent  depths 
derived  from  typical  density  profiles  in  the  tropical  Atlantic  are  60, 
20,  8,  4,  2 cm,  respectively,  for  the  first  five  baroclinic  modes. 

This  places  the  two-laver  svstem  being  considered  about  midway  among 

these . 

3.3.3  Projection  of  the  forcing  function  on  the  vertical  modes 

in  the  previous  subsection,  the  normal  mode  equations,  including 

normal  mode  representations  of  the  forcing  terms,  were  derived.  The 

normtl  mode  forcing  terms  are  now  written  in  terms  of  the  layer  forcing 

terms.  Following  Lighthill  (1969): 

, . N O.H.C^  (F.,  G.) 

(F\  Gk)  = j -T  J 1 -1  1 


o.h  (c1;) “ 

i p j 


w 
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In  particular,  with  F2  = G2  = Fi  = °>  and  Gi  = Gj/O*  a latitudinally 
uniform,  time-dependent,  cross-equatorial  windstress  of  the  form: 


(y) 


Gk  = 


Ck 


PjHj (Ck)2  + 02H2(Ck)2 


Fk  = 0 


For  the  previously  selected  parameter  values,  we  obtain: 


/y) 

80m 


80m 


Thus,  the  forcing  is  the  same  on  each  of  the  two  modes  and  has  a 
magnitude  which  is  equivalent  to  that  when  forcing  a single-layer 
system  80  meters  thick.  For  illustration,  another  case  is  considered. 

A windstress  applied  uniformly  over  both  layers,  with  the  same  densities 
and  thicknesses  as  previously: 


T(y> 

G1  c2  H]  + H2 


results  with: 


G1  = ± 


(y) 


99.4m  ’ 


580  m 
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In  this  case,  most  of  the  forcing  goes  into  the  external  mode,  as  is 
intuitively  apparent. 

3.3.4  Solution  to  the  horizontal  equations 

Once  the  vertical  dependence  has  been  separated  and  harmonic 
dependence  in  longitude  and  time  has  been  introduced,  the  problem  is 
reduced  to  that  of  solving  an  ordinary  differential  equation  in  latitude. 
The  resulting  equation  for  v (u  and  ~ have  been  eliminated)  in  the  homoge- 
nous (unforced)  case  is  (after  Matsuno,  1966): 


? 9 2 ? 

CT  v + (a  - k C + 9C 

g yy  g g 


k 

a 


„ 2 2 

0 y ) v = 0 


where  dependence  of  the  form  e^^ct  l<x'*  ),as  been  assumed  and  = .^gH 
for  the  mode  under  consideration. 

Applying  the  boundary  condition  for  y = + ® where  v -*  0,  the 
equation  admits  to  solutions  of  the  form: 


2C 


v (y)  = cn  e 


* > 

g 


where  H is  the  nth  Hermite  polynomial.  This  solution  requires  that: 
n 


2 , 2 2 a„2  k 

a + k C + 9C  — 

g g c 


BC  (2m  + 1) 
g 
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For  k = 0,  n = 0, 

a.  = /BC  or  T = 2tt/  ,/b'c 

0 •/  g 0 v g 

corresponding  to  the  time  scale  previously  discussed.  Dispersion 
diagrams  for  various  horizontal  modes,  n,  are  presented  in  Fig.  23. 

The  axes  are  labelled  in  dimensional  units  for  each  of  the  vertical 
modes  in  the  chosen  example.  While  this  may  be  more  cumbersome  than 
using  a non-dimensional  dispersion  diagram,  it  is  more  readily  refer- 
enced with  real  quantities. 

The  line  labelled  n = -1  on  the  diagram  represents  the  equatorial 
Kelvin  wave.  This  is  an  additional  solution  to  the  normal  mode  equa- 
tions for  the  case  where  v = 0 everywhere.  Consequently,  the  above 
equation  for  v alone  cannot  provide  this  result.  The  Kelvin  wave  was 
introduced  by  Moore  (1968). 

It  can  be  seen  in  Fig.  23  that  Rossby  wave  modes  have  periods 
in  excess  of  40  days  and,  therefore,  are  not  likely  candidates  for 
explaining  observations  in  the  period  range  from  two  to  three  weeks. 
Since  it  will  later  become  apparent  that  the  Yanai  (n  = 0)  mode  is  of 
greatest  interest,  more  detailed  dispersion  curves  for  these  waves 
are  plotted  in  Fig.  24. 

The  Yanai  wave  is  sometimes  termed  the  mixed  Rossby-gravity 
wave  because,  for  large  values  of  k,  the  westward  propagating  branch 
resembles  a Rossby  wave,  whereas  the  eastward  branch  resembles  a 
gravity  wave.  For  k = 0,  it  lies  in  between.  It  is  possible  that 
the  Yanai  wave  is  in  some  way  related  to  inertial  oscillations  since 
its  period  for  k = 0 is  the  inertial  period  at  y = + Rj . 
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3.3.5  The  Forced  Problem 

The  solution  to  the  general,  linear,  initial-value  problem  will 
not  be  presented  here.  It  is  a very  complex  and  interesting  problem  in 
its  own  right  and  has  been  treated  at  length  by  Cane  and  Sarachick  (1976) 
and  Lighthill  (1969).  It  is  sufficient  to  say  - and  this  point  is 
demonstrated  in  the  numerical  experiment  section  - that  most  of  the 
energy  of  the  response  can  be  accounted  for  by  the  first  several  free, 
horizontal  modes  of  the  system,  plus  a steady-state  component,  when 
present.  The  free  modes  have  been  discussed.  The  steady  solution  is 
quite  simple  to  analyze  since  it  is  the  Sverdrup  transport  resulting 
from  the  east-west  structure  of  the  meridional  forcing. 

3.3.6  Solution  to  the  Steady  Problem 

Moore  (personal  communication)  outlined  the  following  analysis. 
Setting  time  derivatives  to  zero  in  the  linear  equations: 


-Byv1  = -g'H,  -Syv2  = -p2 


I 
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If  k ^ 0,  then  = u^  = v,,  = 0 . 

If  k = 0,  then  v = = 0;  1 | 5 o ; |u^|  ^ 0 . 

For  \ - 2400  km 

2 

r = 1 dyn/cm“ 

H = 40  m 

|vj | = 2.86  cm/ sec. 


3 . 4 Numerical  Experiments 
3.4.1  Numerical  Method 

Taking  RQ=  1 in  Eqs.(3),  Section  3.2  (which  is  equivalent  to 
setting  the  mean  flow  scale,  U,  to  C ) , the  following  equations  result: 


- - 2 3 

u,  = -iku.u,  -v,u,  + yv  - ik + F/0  R, 

1 11111  d 

t y 


2 3 

v = -iku  v -yu,  - + G/P  R. 

1 1111  d 

t v 


\ \ = -ikui  (YV  V -vi(Y'i  - 2 

t t y y 


* ) - (1  + y\-  \)  (iku.+v 


rvi  ) 
y 


V 
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-iku2u2  - v 5 + yv9  - ikC^  + 1^) 

y 


-ikuv9  - yu2  - ^ - il2 

y y 


5"^  + a\  = iku,  (5^  4-  c/^)  - V9  (5^  + a-T^  ) 

t t y y 


(H  + 5^  + orT^)  (iku9  + vo  ) 

“y 


Notice  that  the  third  and  sixth  equations  are: 


v\ 


+ ^ « B 


which  results  in: 


T!j  = (B  + oA)  / (5  + cry) 

\t  = (yB  - JA)  / (6  + cry) 

Conditions  at  zonal  boundaries  were  treated  as  follows: 


ult  = -ik\ 

U9  = - ik  (T^  + 

± yBul  - 

1 

f=^ 

1 

II 

CM 

PQ 

+ 1 

y 

y 

vi  = v2  = o 

at  y = ± yB 
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This  results  in  a system  of  six  equations  of  the  form: 

^ (*")  , * = 1,6 

where  F involves  9 and  its  derivatives  in  y. 

It  was  found  that  a fourth-order  differencing  scheme  for  deriva- 
tives in  y was  both  efficient  and  accurate  (Cane  1975).  Time  integra- 
tion was  done  by  the  "leapfrog-trapezoidal"  method  (Haltiner , 1971) . 

These  techniques  are  now  described. 

Space  differencing: 

f “ C8<51+1  - »i-l>  - «l+2  - Si-2>’  1 

In  cases  where  i is  within  1 grid-point  of  the  boundary,  a slightly 
different  form  is  necessary: 

~ « (6«3  - 3«2  - - 2?1)  / 6 Ay 

39 

w (-993  + 1R?2  + 2$^  - 119^  / 6 Ay 

Boundary  conditions  were  managed  as  follows: 

(Mean  flows  were  assumed  to  be  zero  at  y = yb . Tendencies  for  u^,  u2 
were  computed.  Following  time  integration,  T^,  were  determined  on 
the  boundary  in  terms  of  adjacent  values  °f  and  the  new  values 

of  u^  and  u2  on  the  boundary. 

The  "leapfrog-trapezoidal"  method  is  a two-step  or  predictor- 
corrector  method.  It  was  chosen  for  efficiency  and  because  it  does  not 
adversely  affect  phases  and  amplitudes  of  wave-type  solutions.  For  long 
integrations  or  those  involving  frictional  processes,  however,  other 
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schemes  may  be  preferable.  The  method  is  summarized  as  follows: 


let 


r ( $k> 


where  k represents  the  time-step  and  F the  tendency  function. 
We  define  a predicted  value: 

•A. 

? = «k_1  + 2Atr(5k) 


where  fat  is  the  time-step  size. 
The  corrected  or  final  value  is: 


.k+1 


= 


— r r 
2 L 1 


■k 

($  ) 


ri5k)  ] 


This  scheme  also  has  the  desirable  characteristic  of  damping  the  com- 
putational mode. 

After  some  experimentation  a grid  spacing  in  latitude  of  28  km 
was  chosen.  Linear  stability  of  the  calculation  required  that  the  time- 
step  be  less  than  about  3 hours,  so  a value  of  2 hours  was  chosen  and 
found  to  be  satisfactory.  With  these  step  sizes  an  integration  of  2 
months  (model  time)  required  only  about  15  minutes  of  CPU  time  and 
very  little  core  on  a UNIVAC  1106. 

The  six  dependent  variables  computed  were  decimated  to  four 
steps  per  model  day  before  storage  on  disk  and  tape  for  subsequent 
plotting  and  analysis. 

The  forcing  function  employed  was  essentially  uniform  in  lati- 
tude. It  was  tapered  off  inside  the  zonal  boundaries  to  prevent  the 
immediate  generation  of  gravity  waves  which  might  contaminate  the 
equatorial  solution.  Time  dependence  was  similar  to  the  latitude 
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dependence,  being  a hyperbolic  tangent  function.  The  form  of  the 
forcing  was: 

G = T/Ltanh  a(t-t^)  - tanh  aCt-t^)]  * [tanh  b(y-y^)  - tanh  b(y+y^)]/  4 

Where  y^  is  the  extent  in  latitude,  north  and  south,  and  t^  and  t^  are 
the  onset  and  decay  times,  respectively.  r'  is  the  scaled  forcing 
amplitude.  This  form  allows  for  the  application  of  a pulse-type 
forcing  of  any  duration.  For  example,  for  t^  60  days,  near  zero  and 
a sufficiently  large,  G approaches  a step  function  in  time.  The  form  of 
the  forcing  will  be  shown  (Fig.  25)  in  the  individual  examples. 

The  mean  flow  was  taken  to  be  a Gaussian,  zonal  flow  in  each 
layer  specified  by  the  parameters  and  A . 

_ 9 9 

u^  = Up  exp(-y  /2A.)  i = 1,  2 

i 

^ (the  pressure  field)  was  chosen  to  be  in  geostrophic  balance  with 
u (recall  the  third  and  seventh  equations  in  (2))  . Layer  thicknesses 
were  specified  at  the  equator  so  that,  in  the  presence  of  mean  flow, 

T]  increases  from  zero  at  the  equator  (see  Fig. 26).  were  chosen 

to  be  consistent  with  the  observed  latitudinal  extents  of  the  equatorial 
currents,  thus  mean  flows  were,  subsequently,  assyptotically  zero  well 
within  the  zonal  boundaries. 

2 lot 

Since  longitudinal  dependence  of  the  form  e has  been  assumed, 
all  the  dependent  variables  are  complex.  This  was  achieved  by  using 
complex  FORTRAN  variables.  Real  and  imaginary  parts  of  the  model  output 
were  combined  appropriately  in  later  analysis  to  recover  features  of  the 
longitudinal  structure  such  as  zonal  propagation.  The  forcing  function 
.s  assumed  to  have  only  a real  part  which  is  equivalent  to  a cos(kx) 


I 
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dependence.  The  forcing  is  thus  a stationary  or  standing  wave,  depen- 
ding on  the  temporal  dependence. 

3.4.2  Cases  Studied 

The  numerical  experiments  consist  of  six  cases.  The  density 
structure  in  all  cases  is  as  described  in  Chapter  3,  Section  1.  The 


are 

as  follows: 

1. 

K = 0 

Forcing  = Step  1 

2. 

K = 2 rr/2400 

km 

Forcing  = Step  1 

3. 

K = 2 TT  / 2400 

km 

Forcing  = Step  2 

Zero 

mean  current 

4. 

K = 2 n / 2400 

km 

Forcing  = Pulse 

5. 

K = 2 n / 2400 

km 

Forcing  = Step  2 

Mean 

current  present 

6. 

K = 2 rr  /2400 

km 

Forcing  = Pulse 

Figure  25  shows  the  time-dependence  for  the  three  forcing  func- 
tions used.  The  amplitude  or  maximum  value  of  windstress  is  chosen  to 
2 

be  0.25  dyn/cm  . Since  this  is  a linearized  model,  i.e.,  the  non- 
linear terms  are  simulated  by  variable  coefficients,  the  amplitude  of 
the  response  is  directly  proportional  to  the  amplitude  of  the  forcing. 
The  response  depends  in  no  other  way  on  the  forcing  magnitude. 

Characteristics  of  mean  flow  are  given  below  for  cases  where  it 
is  indicated. 

Uq  = 40  cm/sec.  = 300  km 

Uq  = 40  cm/sec.  ^2  = km 


Horizontal  profiles  of  these  flows  are  shown  in  Fig. 26.  In  all  cases. 
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integration  was  carried  out  for  sixty  days  of  model  time. 


3.4.3  Derived  Quantities 

In  addition  to  the  six  dependent  variables,  several  other  pertinent 
quantities  were  calculated.  These  include  the  internal  and  external 
modes  of  the  v-component  for  case  4,  total  u-component  for  cases  5 and  6, 
and  near-equatorial  v-components  as  a function  of  longitude  and  time. 

The  derivations  are  now  described. 

The  external  and  internal  modes  were  found  from  the  appropriately 
weighted  sum  of  v-components  in  the  two  layers.  From  Section  3.3,  we 
have: 


v*  (ext)  - (v 


-2 


! ' 2 


v (int)  = (v ^ - ,|2v2)/  2 
The  total  u-components  were  found  by  adding  the  fluctuation  u- 
components  to  the  prescribed  mean  flows  ; 

% = “i + v ni  - v 

% = G2  + U2r ’ G2  - U2. 

where  u^ , u^  are  real  and  imaginary  parts  of  u respectively. 

v-components  as  functions  of  longitude  and  time  were  derived 
using  averages  over  a finite  latitude  band: 


V = - 2NTT 


N 

r* 


V 


r 


J Vy 


(i  - 1,  2) 


where  y = y = 280  km  — 2.5°, 
° N 
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The  resulting  time  series  of  ( ) were  then  multiplied  by  the 
harmonic  longitude  dependence,  and  the  real  part  taken: 

vi(x,t)  = (v^  )r  cos(kx)  - ( sin(kx)  (i=l,2) 

where  k = 2tt/2400  km;  x = x^  , j = 1,121 

X121_X1  = km  or  two  cycles. 

and  dx  = x,.,-x.  = 40  km 
J+l  J 

3.4.4  Presentation  of  Results 

Results  from  cases  3 and  6 are  shown  in  Figs.  27-64.  Cases  1 
and  2 are  discussed  but  plots  of  output  are  not  included;  case  2 is 
not  qualitatively  different  from  case  3;  and  case  1 is  a test  run  with 
k=0.  Most  of  the  indicated  plots  are  latitude-time  maps  of  real  and 
imaginary  parts  of  the  six  dependent  variables,  the  total  u-component 
and  the  vertical  modes.  The  remainder  (Figs.  33,  34,  43,  44,  53,  54, 
63  and  64)  are  longitude-time  representations  of  near  equatorial  v- 
components . 

All  variables  are  presented  for  completeness.  Most  discussion, 
however,  concerns  the  meridional  components.  Latitude-time  maps 
extend  from  5°S  to  5°N  and  cover  the  entire  sixty-day  integration 
time.  Figure  65  is  a second  plot  of  the  real  part  of  v^ , case  4,  but 
covering  the  entire  latitude  domain.  With  the  exception  of  some 
artificial  boundary-generated  noise  near  the  end  of  the  integration, 
the  response  is  limited  latitudinally  to  about  +5°.  The  larger 
radius  of  deformation  for  the  system  is  267  km  or  about  2.4°  which  is 
indeed  consistent  with  the  trapping  scale  of  the  model  results.  There 
are,  in  fact,  extra-equatorial  motions  present  in  the  response,  e.g., 
inertial  oscillations,  but  these  are  at  least  an  order  of  magnitude 


I 


6? 

smaller  in  amplitude  t!  ■:  the  equatorial  response  For  these  i . : : , .^n. 
the  +5°  scale  was  chost.n  s : nr  • it  provides  • more  J tailed  presentation 
In  all  the  Jatitudu-t Jr  • plots  one  can  n-.e  pn  . . r ed,  si  vu  • : - 
week  oscillations.  n the  i'  t.gitudi  -time  1 t.  , .i  J inir.  ml , westward 

propagation  is  evident.  A detailed  u in  tuo.,  i •>..  > ' t . : > .on 

in  Chapter  4. 
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4.  DISCUSSION  OF  RESULTS 

4 . 1 Comparison  of  Model  Results  with  Linear  Theory 

To  ascertain  that  the  numerical  scheme  made  sense,  four  runs  were 
made  with  zero  mean  flow  (Cases  1-4,  Section  3.4.2) „ The  results 
of  these  runs  were  projected  onto  the  theoretical  free  modes  of  the 
system  to  determine  if  the  model  was  consistent  with  linear  theory. 

For  Case  1,  k = 0;  for  Cases  2 to  4,  k = 2t/2400  km,  The  resul- 
ting v-components  were  first  combined  to  form  vertical,  normal  modes, 
and  the  latter  then  projected  onto  the  horizontal  modes,  corresponding 
to  the  n = 0,  2,  4,  6,  8,  10,  12,  and  14  Hermite  functions.  Projec- 
tions on  odd-order  modes  are  zero,  as  is  to  be  expected  for  the 
prescribed  forcing.  In  Fig.  67,  the  percentage  of  variance  accounted 
for  by  each  of  these  modes  is  presented.  These  results  show  that  the 
majority  of  the  response  is  confined  to  the  n = 0 mode,  or  Yanai  wave, 
with  percentages  tapering  off  rapidly  as  n is  increased.  The  Hermite 
function  coefficients  which  result  from  this  analysis  are  time  series, 
and  have  been  further  analyzed  to  extract  periods  and  phases  in  time. 
The  results  are  in  excellent  agreement  with  the  dispersion  relations 
derived  in  Chapter  4. 

The  analysis  procedure  is  now  described. 

k <t> 

Let  v = — 2rr-  = Sverdrup  transport, 
s °l  1 

T 

max 

1 f* 

where  (t>  = - — v dt  (i.e.,  the  average  forcing). 

T * « 

max  0 


w 
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Then  for  k 4 0 , 

vx  (y,  t)  s vl  (y , t)  + i (v  (y,C)  + vg)  ® v1  (y,c)  + iv  (y,t) 
r i r l 

v2(y,t)  = v2  (y , t)  + iv2  (y,t)  ; 
r 1 

consequently, 

i • ii 


as  was  found  earlier. 


Forming  the  vertical,  normal  modes: 

~ 1 „ 1 / , „ 1 . , 
v = v2  + C2  V1  external 

-2  2 2 

v = C ^ v'  + C2  v2  internal  , 

where  the  C's  represent  the  eigenvectors  of  the  normal  mode  problem. 

Each  of  these  is  then  projected  on  the  appropriate  Hermite 
functions : 


00 


(t) 


v dy 


Y£  = e H < -$-) 

R , 


t = 1,2 


y = (-«,  ”) 


is  such  that  Y is  assymptotical ly  zero  at  y^  = + 15°,  so  the  limits 
of  integration  are  placed  here  without  much  loss  of  accuracy.  The 

l 

C (x)  constitude  time  series  of  Hermite  function  coefficients, 
m 


L 
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The  variances  in  Fig.  26  were  comouted  as  follows: 


l 

D 

m 


K <[»2 

( v (y , t) ) 


dy  dt 


For  k = 0,  these  time  series  are  sinusoidal.  The  periods  agree 
with  those  predicted  by  the  dispersion  curves  unequivocally  since, 
for  each  horizontal  (and  vertical)  mode,  there  is  only  one  frequency. 

For  k 4 0,  however,  the  series  are  not  simple  sinusoids,  but  rather, 
composite,  periodic  signals.  Die  imaginary  parts  of  these  series 
(which  were  zero  everywhere  for  k = 0)  contain  components  of  the  same 
frequencies  but  of  differing  phase  from  those  in  the  real  parts.  It 
was  necessary,  therefore,  to  extract  the  frequencies  and  phases  of 
these  components  in  order  to  determine  the  longitudinal  propagation 
and  dispersion  characteristics. 

The  dispersion  relation  (Fig.  23)  implies  that,  for  k i-  C , there 
should  be  three  frequency  components  for  n ^ 2,  two  for  n = 0,  and  one, 
the  Kelvin  wave,  for  n = -1;  i.e.,  for  n 2:  0 , there  is  an  eastward  and 
a westward  gravity  wave,  as  well  as  a Rossby  wave.  As  stated  above, 
we  will  not  consider  the  Kelvin  wave.  The  numerical  solution,  after 
integrating  out  the  structure  in  latitude  as  described,  may  be  represen- 
ted by: 

v(x,t)  = Re  | (v^  (t)  + iVj  (t))  e f 

where  v^  and  v^  are  the  real  and  imaginary  parts  of  the  model  results 
projected  on  a mode.  Assuming  that  v^  and  v^  are  linear  combinations 


I 


of  harmonic  functions  in  time  with  frequencies  cr^  and  cr^ : 
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112  2 

v = A cos(at)  + B sin(cr.t)  + A cos(a0t)  + B sin(a„t) 
r r 1 r 1 r 2 r 2 

112  2 

v^  = A^  cos(cjjt)  + B^  sin^t)  + A^  cos(a,,t)  + B^  sin(a2t) 

Thus,  for  each  frequency,  a.,  we  have,  after  some  manipulation: 

v(x,t).  = I (A^  - B-))  cos(a.t  - kx)  + (B^  + A^)  sin (a.t  - kx)  + 

' ’ j 2 L r i 3 r i j 

(A^  + B^)  cos  (d.t  + kx)  + (B''  - A^)  sin  (a.t  + kx)  j 
ri  3 ri  j J 

The  first  two  terms  correspond  to  eastward  propagation  and  the  last 
two  to  westward  propagation.  The  coefficients  A^ , B^,  etc.  were 
evaluated  from  the  model  results  by  the  method  of  least  squares. 

For  n = 0,  there  are  only  two  frequencies  present.  For  n s 0, 
however,  we  also  have  a Rossby  wave.  Only  the  two  inertio-gravity 
wave  frequencies  were  explicitly  extracted  since  the  Rossby  waves  were 
of  periods  as  large  as  and  larger  than  the  integration  time.  Subse- 
quent visual  comparison  of  original  and  reconstituted  time  series 
ascertained  that  virtually  all  short-period  variance  had  been  accounted 
for  by  and  , and  that  the  remainder  was  a long-period  variation 
consistent  with  the  Rossby  wave. 

The  numerical  solution  (no  mean  flow)  thus  consists  of  a combi- 
nation of  the  free  modes  of  the  system  with  the  correct  frequencies 
and  propagation  directions.  Moreover,  the  greater  part  of  the  variance 
of  the  response  is  contained  in  the  lowest  order  horizontal  modes 
(Fig.  67).  A typical  response  amplitude  is  20  cm/sec  for  v. 

Examination  of  Fig.  67,  in  conjunction  with  the  three  forcing 
functions  (Fig.  25),  reveals  that  the  response  to  the  pulsed  wind  is 
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confined  almost  entirely  to  n = 0 for  two  modes.  For  both  stepped- 
wind  cases,  there  is  a broader  distribution  of  variance  over  horizon- 
tal modes.  In  particular,  eight  horizontal  modes  account  for  only  727. 
of  the  variance  of  the  internal  mode.  This  broadening  can  be  attributed 
to  the  existence  of  many  Rossby  wave  modes  to  satisfy  initial  conditions 
in  the  presence  of  the  Sverdrup  transport.  Only  in  the  stepped-wind 
cases  is  this  transport  non-negl igible . For  Case  1,  where  k = 0,  707, 
of  the  variance  is  found  in  the  internal  mode.  This  is  consistent  with 
the  other  stepped-wind  cases  in  that  the  internal  mode  is  excited  more 
than  the  external  mode  by  the  steady  component  of  the  wind.  Internal 
mode  frequencies  are  lower  than,  and  thus,  closer  to  the  energy  peak  in 
the  steady  wind.  This  effect  is  in  addition  to  the  steady-flow  effect, 
which  does  not  exist  in  Case  1.  The  reader  is  reminded  at  this  point 
that  the  prescribed  forcing  is  projected  equally  on  the  baroclinic  modes. 

Case  4 is,  as  a result  of  these  considerations,  the  preferred  case 
for  study  of  the  short-term  response  because  virtually  all  the  variance 
is  confined  to  two  lowest  order  horizontal  modes,  and  the  intermittent, 


pulse-like 

forcing 

is  probably  more 

representative  of 

real  conditions. 

Table  5 

Yanai  Wave  Characteristics  for  Model  Modes 

Period 

Phase  Speed 

Amplitude 

(d imen  s ion 1 e s s ) 

External 

westward 

(days) 

(cm/sec) 

Case  2 

Case  3 

Case  4 

Case  1 

16.7 

166 

0.277 

0.189 

0.294 

eastward 

8.4 

329 

0.063 

0.039 

0.072 

Internal 

westward 

23.1 

120 

0.177 

0.244 

0.297 

eastward 

14.7 

188 

0.073 

0.094 

0.158 

External 

11.9 

0.103 

Internal 

18.5 

0.163 
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Table  5 contains  amplitudes  of  Yanai  wave  components  for  each  of 
the  four  cases  discussed.  For  Case  1 (k  = 0) , these  amplitudes 
merely  reflect  the  previous  argument  concerning  internal  and  external 
modes  since  there  is  only  one  wave  in  each.  In  the  other  cases,  how- 
ever, there  is,  in  addition,  a systematic  difference  between  the  ampli- 
tudes of  westward-  and  eastward-propagating  waves:  the  amplitudes  of 

the  westward  waves  are  always  the  larger.  This  can  be  explained  by 
inspection  of  the  dispersion  diagrams  (Fig.  23)  and  Table  5 and  by  noting 
that  the  phase  speed,  a/k,  is  about  twice  as  large  for  the  eastward 
branch  of  the  Yanai  wave.  Consequently,  the  westward  wave,  with  slower 
phase  speed,  tends  to  absorb  more  energy  before  it  propagates  suffi- 
ciently far  to  be  out  of  phase  with  the  stationary  forcing. 

In  summary,  of  all  components  appearing  in  the  response,  the 
dominant  one  is  the  Yanai  wave,  particularly  the  westward-propagating 
branch,  with  both  external  and  internal  modes  of  similar  magnitude. 

The  longitude-time  plots  for  Cases  3 and  4 (Figs.  33,  34,  43,  and  44) 
show  a dominant,  westward-propagating  disturbance  superimposed  by  beat 
patterns  related  to  the  frequency  differences  between  internal  and 
external  modes  and  by  interference  from  eastward-propagating  waves.  In 
Case  3 (Fig.  33),  in  the  upper  layer,  the  picture  is  further  complicated 
by  the  presence  of  steady  flow  (Sverdrup  flow)  and  Rossby  waves. 

Phase  speeds  were  extracted  for  the  dominant,  westward-  and  east- 
ward-propagating features  directly  from  the  plots  (Table  6).  These 
values  agreed  with  those  of  the  dominant,  theoretical  mode  in  each  case. 

Cases  5 and  6,  which  include  mean  current,  are  not  analyzed  by  the 
foregoing  techniques  because  vertical,  normal  modes  cannot  be  defined 
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in  the  same  way  . It  is  necessary,  therefore,  to  make  direct  compari- 
sons of  plots  of  output  to  estimate  the  effects  of  the  advective  terms. 


Table  6:  Phase  Speeds  (cm/sec)  Extracted  from 

Longitude-Time  flaps;  Pulsed  Forcing 


Case  4 Case  6 

(no  mean  flow)  (mean  flow) 


Upper 

Lower 

Average 

Upper 

Lower 

Average 

Westward 

152 

142 

147 

186 

159 

173 

Eastward 

200 

186 

193 

174 

186 

180 

4 . 2 Effects  of  Mean  Current 

Examination  of  longitude-time  plots  for  Cases  5 and  6 (Figs.  53, 
54,  63,  and  64)  indicates  that  westward-propagating  disturbances 
dominate  even  more  than  in  linear  cases.  Average  westward  phase 
speeds  tend  to  be  increased,  while  eastward  phase  speeds  tend  to  be 
decreased.  This  difference  may  constitute  a Doppler  shift  associated 
with  the  current  in  the  upper  layer.  The  reason  that  the  Undercurrent 
has  a much  smaller  effect  in  this  shift  is  that  it  is  narrower  than 
the  smaller  of  the  two  radii  of  deformation,  and  thus,  has  only  a 
secondary  Doppler-shifting  effect,  while  the  width  of  the  current  in 
the  upper  layer  is  greater  than  the  larger  R^.  Phase  speeds  derived 
from  x-t  plots  for  Cases  3 to  6 are  recorded  in  Table  6.  The  Doppler- 
shift  described  above  tends  to  cause  the  frequencies  of  the  eastward 
and  westward  waves  to  be  closer  together.  These  results  should  be 
viewed  as  qualitative  since  the  average  phase  speeds  derived  from  the 
plots  are  only  approximations.  In  short,  the  Doppler  effect  produces 
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a signal  which  is  closer  to  containing  a single  frequency.  Indeed, 
in  the  latitude-time  maps  of  lower-layer  v-com ponents  (Figs.  32,  39, 

50,  and  60),  the  signal  seems  to  be  more  nearly  simple  harmonic  in 
time  than  in  the  linear  cases.  However,  while  some  apparent  Doppler 
effects  have  been  noticed,  the  concept  is  limited  in  this  situation 
since  the  latitudinal  scale  of  the  waves  is  of  the  same  order  as  that 
of  the  zonal  currents. 

The  amplitudes  in  the  mean  flow  cases  tend  to  be  greater  than  in 
the  other  cases.  This  implies  an  initial  transfer  of  momentum  and 
energy  from  the  mean  current  to  the  fluctuations,  but  does  not  neces- 
sarily constitute  an  instability.  An  instability  would  require  that  the 
amplitudes  grow  with  time.  It  might  be  said  that  the  interaction  with 
the  mean  current  is  an  amplification  of  the  forcing  signal.  It  is 
likely  that  the  proper  choice  of  mean  flow  parameters  would  yield  an 
instability  for  the  forcing  and  wavelength  used.  This  problem  has  been 
addressed  by  Philander  (1976)  for  a very  similar  system,  but  will  not 
be  examined  here. 

The  immediate  effects  of  the  mean  current  through  non-linear  inter- 
action can  best  be  seen  in  Figs.  29,  36,  47,  and  57.  The  signs  of  vu^ 
are  such  that  negative  fluctuations  in  u (north  of  the  equator)  in  the 
upper  layer  and  positive  fluctuations  in  the  lower  layer  are  enhanced, 
while  fluctuations  of  the  opposite  signs  are  suppressed. 

The  presence  of  mean  zonal  currents  with  horizontal  shear  modify 
the  responses  in  the  cases  considered,  but  the  picture  is  still  one  of 
equatorially  trapped  waves  with  much  the  same  structure  as  in  cases 


\ 


71 


without  mean  flow.  These  results  suggest  that  the  system  responds 
with  periods  and  spatial  scales  which  are  determined  primarily  by  the 
basic  system  parameters,  namely,  the  density  structure.  While  it  is 
likely  that  values  of  east-west  wavenumber  and  mean  flow  parameters 
can  be  found  for  which  the  response  is  entirely  different,  e.g., 
unstable,  small  deviations  from  the  chosen  values  will  have  much  the 
same  results  as  those  presented. 

4 . 3 Comparison  of  Observations  with  Theoretical  and  Numerical  Results 

As  discussed  in  Chapter  2,  considerable  caution  is  required  in  the 
inference  of  wave  parameters  from  GATE  observations.  The  implications 
are  actually  that  certain  scales  of  motion  are  observed  more  frequently 
than  others.  Even  the  results  of  one  of  the  previously  discussed 
model  runs,  vith  cr  without  mean  flow  if  sampled  in  the  same  fashion 
as  GATE  observations  were  made,  would  provide  somewhat  ambiguous 
results.  With  these  limitations  in  mind,  the  following  arguments  are 
presented . 

Since  the  best  estimate  of  a wave  parameter  from  the  observations 
is  the  16-day  period,  marks  on  the  a-axes  (see  Fig.  23)  corresponding 
to  this  value  show  that,  in  the  external  mode,  only  the  westward- 
propagating  Yanai  wave  and  the  Kelvin  wave  curves  are  intersected. 

For  baroclinic  modes  whose  equivalent  depths  exceed  that  of  the  exter- 
nal mode  of  the  model  (H  = 27.3  cm)  up  to  a maximum  that  would  be 
expected  in  the  tropical  ocean,  say  about  200  cm,  this  result  holds. 

, for  the  lower  baroclinic  modes,  a 16-day  wave  would  have  a 
i ..  peed  between  about  120  and  200  cm/sec  and  a wavelength  between 
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1700  and  2700  kin.  For  higher  baroclinic  modes,  such  as  the  internal 
mode  of  the  model,  the  eastward-propagating  Yanai  wave  becomes  a 
possibility.  ITiese  results  are  summarized  in  Fig.  24,  which  contains 
plots  of  period  vs.  phase  speed  for  several  equivalent  depths  for  the 
Yanai  wave.  Also  plotted  are  curves  of  constant  wavelength.  The 
rectangle  superimposed  on  the  curves  for  the  westward  wave  represents 
the  ranges  of  these  quantities  inferred  from  the  observations  for  the 
n = 0 case  in  Table  5.  The  plot  suggests  that  the  observed  wave  para- 
meters are  consistent  with  the  first  or  second  baroclinic  mode  in  the 
continuously  stratified  case  and  for  the  external  mode  of  the  2 ' - ! aver 
system  under  consideration.  Phase  speeds  for  other  values  of  n in 
Table  5 are  too  low  to  be  in  agreement  with  linear  theory.  It  is 
possible,  however,  that  an  eastward-propagating  Yanai  wave  of  smaller 
amplitude  is  present  and  may  be  responsible  for  the  assymmetry  observed 
in  the  y-t  maps  (e.g..  Fig.  3).  Indeed,  when  the  results  of  the 
numerical  integrations  were  analyzed,  such  was  found  to  be  the  case. 

As  shown  in  Section  2.2,  there  are  indications  that  shorter  period 
fluctuations  exist  in  the  observations.  These  appear  to  have  periods 
between  5 and  10  days.  Figure  23  shows  that  such  periods  are  consistent 
with  n = 2 iner tio-gravity  waves  for  a wavelength  of  from  2000  to 
2500  km  for  the  first  several  baroclinic  modes  and  the  external  model 
mode.  Consideration  of  the  Kelvin  wave,  at  least  for  the  major  fluctua- 
tion, is  ruled  out  since  it  requires  v to  be  zero.  In  fact,  all  odd- 
order  Hermite  functions  are  zero  at  the  equator.  Unfortunately,  phase 
speeds  cannot  be  estimated  for  these  shorter  period  fluctuations  from 
the  observations.  Sea- level  measurements  in  the  equatorial  Pacific 
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analyzed  by  Wunsch  and  Gill  (1976)  show  the  existence  of  first  baro- 
clinic  mode  equatorial  waves  with  periods  of  the  order  of  four  days. 

Wunsch  and  Gill  also  present  an  argument  for  the  existence  of  meri- 
dionally forced  Yanai  waves. 

In  summary,  it  has  been  shown  that  the  wave  parameters  estimated 
from  the  observed  current  records  are  consistent  with  the  free  solu- 
tions of  the  linear,  equatorial  wave  problem  for  at  least  the  westward- 
propagating  Yanai  wave  for  low-order  baroclinic  modes. 

Comparison  of  model  results  with  and  without  mean  current  indicates 
that  the  response  is  wavelike  in  either  case,  thus  making  the  comparison 
with  linear,  free  modes  more  meaningful.  Comparison  of  observations 
with  model  results  where  the  mean  flow  is  included  is  also  necessary, 
however.  One  case  is  selected  for  comparison  as  being  most  representa- 
tive of  real  conditions. 

Previous  descriptions  (Section  2.4)  of  observed  winds  during  GATE 
imply  that  the  pulsed  wind  (Cases  4 and  6)  is  probably  the  simplest  realistic 
approximation  to  real  conditions  for  examining  short-term  responses.  In 
any  case,  it  has  been  shown  that  a wavelike  response  results  for  step- 
like forcing  as  well  as  the  pulse. 

The  model  results,  in  conjunction  with  the  theoretical  analysis 
of  the  linear  problem,  suggest  that  the  oscillations  of  the  Undercurrent 
constitute  a north-south  advection  by  a predominantly  wavelike  process. 
This  conclusion  would  be  apparent  in  the  observations  if  the  north- 
south  or  trapping  scale  of  the  v-componont  were  much  greater  than  the 
width  of  the  Undercurrent.  However,  based  on  observations  alone,  it 
is  not  obvious  whether  the  oscillations  are  primarily  wavelike,  as 
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described,  or  are  manifestations  of  an  inertial  jet  phenomenon.  'Ihe 
latter  hypothesis  is  weakened  if  the  zonal  flow  is  assumed  to  be 
geostrophically  balanced,  and  the  pressure  field,  as  well  as  the  flow 
field,  experiences  the  same  initial  disturbance;  i.e.,  the  lateral, 
Lagrangian  motion  of  a free  jet  might  be  described  by  the  following: 

dv 

c . £. 


where  u^  is  the  eastward  speed  of  the  jet  and  v is  the  northward 
velocity  of  the  position  of  the  core. 

For  u^  = 80  cm/' sec. 


T = 


= 17  days. 


If,  however,  u„  is  the  geostrophic  part  of  u„ , then 
Og  ° 11  0 


8yu0g  = -py  * 


which  greatly  increases  T if  u - u is  taken  to  be  a perturbation. 

0 Og 

This  argument  is  presented  as  a caution  against  an  over-simplified, 
free-jet  hypothesis. 

A point  concerning  the  actual  forcing  is  worth  mentioning.  Ln 
the  model,  the  system  is  initially  at  rest,  so  all  subsequent  motion 
can  be  attributed  to  the  known  applied  forcing.  In  the  observations, 
however,  a small  segment  of  an  ongoing  process  is  being  examined,  and 
our  attempts  to  correlate  the  local  winds  with  the  currents  may  be 
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somewhat  naive. 


If  the  observed  fluctuations  are  indeed  produced  by 


the  mechanisms  implied  by  the  model,  it  is  more  likely  that  they  were 
excited  by  a pulse  or  periodic  wind  further  to  the  west  and  previous 
to  the  current  observations.  This  is  implicit  in  the  model  results 
since  the  most  energetic  part  of  the  lower  layer  response  is  initially 
in  the  imaginary  part,  or  ^-wavelength  away  from  the  peak  in  the 
forcing.  Thus,  at  any  given  time  of  observation,  the  response  at  a 


given  location  is  probably  a combination  of  the  effects  of  local 
forcing  and  previously  excited  waves  which  propagated  into  the  area. 

A comparison  of  v-components  (Figs.  13,  14,  59,  and  60)  shows 
that  amplitudes  and  times  of  zero-crossings  are  in  remarkable  agree- 
ment. It  may  be  re-emphasized  that  the  choice  of  parameters,  both  for 
the  density  structure  and  for  the  mean  current,  is  not  critical  to  this 
agreement.  The  choices  were  based  on  observations  and  adjustments 
were  not  made.  The  total  u-cotnponents  (Figs.  13,  14,  61,  and  62)  are 
also  very  similar.  In  particular,  the  maximum  displacement  of  the 
Undercurrent  core  is  seen  to  be  nearly  the  same  in  both  observations  and 
model  results.  In  the  latter  case,  this  depends  to  some  extent  on  the 
magnitude  of  the  forcing,  but  is  affected  by  P.^  as  well.  Due  to  the 
fact  that  the  mean  fields  in  the  model  are  fixed,  no  changes  in  u or  T) 
result  at  the  equator  for  the  prescribed  forcing.  Thus,  shifts  of  the 
Undercurrent  are  represented  by  changes  in  u off  the  equator. 

Figures  66  a and  b shew  equatorial,  meridional  transports  calcula- 
ted from  the  imaginary  parts  of  the  v-components  in  the  model  output, 
with  and  without  mean  current.  These  are  compared  with  plots  of  the 
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same  calculation  done  for  the  observations  (Fig.  191.  This  compari- 
son shows  the  similarity  in  the  phase  relations  of  the  two  layers, 
which,  in  turn,  suggests  similarity  in  the  dynamic  processes. 

In  general,  amplitudes  in  the  observations  tend  to  be  higher 
than  those  in  the  model.  This  is  likely  to  be  due  to  the  choice  of 
windstress;  thus,  the  discrepancy  need  not  be  of  undue  concern.  Hie 
difference  is  more  pronounced  when  comparing  vertical  motions.  An 
estimate  of  vertical  motion  implied  by  the  model  results  can  be  found  by 
examining  Fig.  52.  The  pattern  is  periodic,  just  as  with  the  observa- 
tions, and  the  magnitude  is  simply  the  time-derivative  of  the  interface 
elevations.  (Advective  effects  cause,  at  most,  a factor  of  two 

difference.)  Where  the  observations  imply  vertical  motions  of  the 

-2  -3 

order  of  5 x 10  cm/sec,  the  model  implies  values  nearer  to  10  cm/sec. 

This  may  be  partially  due  to  the  oversimplification  of  the  shear  zone 
by  the  two-layer  approximation.  Another,  more  likely  explanation  is 
that  there  may  be  additional  waves  present  in  the  observations,  e.g., 
Kelvin  waves,  which  have  a maximum  in  w at  the  equator.  For  even  modes, 
w is  zero  at  the  equator.  For  example,  in  cases  where  Kelvin  and  Yanai 
waves  are  additive,  considerably  higher  magnitudes  in  w would  result. 

4 . 4 Future  Research  Suggested  by  the  Present  Investigation 

During  the  course  of  this  investigation,  primarily  in  the  latter 
stages  when  it  became  clear  that  the  response  remained  wavelike  in  the 
presence  of  a mean  current,  certain  new  questions  arose.  Some  of  them, 
such  as  that  concerning  longitudinal  propagation  and  the  effects  of 
zonal  boundaries,  cannot  be  addressed  with  the  model  in  its  present 
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form;  i.e.  , explicit  x-dependence  would  have  to  be  introduced.  On 
the  other  hand,  there  are  a few  questions  which  are  natural  extensions 
of  this  work  which  are  now  briefly  discussed. 

One  question  which  arose  during  the  GATE  Miami  Workshop  (March, 
1977)  concerned  dispersion  characteristics  in  the  presence  of  mean  flow. 
Although  for  certain  mean  flow  profiles  a perturbation  theory  might  be 
developed  to  analytically  derive  a dispersion  relation,  a more  direct, 
empirical  approach,  which  could  be  made  with  the  model  in  this  study, 
is  to  run  the  model  for  a series  of  values  of  the  east-west  wavenumber 
and,  subsequently,  to  determine  frequencies  by  spectral  analysis  of 
the  results.  This  would  be  done  for  several  mean  flow  configurations 
and,  at  first,  an  impulse-like  forcing  function  to  suppress  longer- 
period  Rossby  waves. 

A second  problem  of  interest  is  that  of  the  energetics  of  the 
system.  Cane  (1975)  addressed  this  problem  in  his  model  calculations, 
providing  an  indication  of  how  it  should  be  approached.  It  would  be 
useful  to  calculate  the  total  (spatially  integrated)  accumulated 
kinetic  and  available  potential  energies  as  a function  of  time. 
Juxtaposed  with  these  should  be  calculations  of  cumulative  work  done, 
as  a function  of  time,  by  the  various  flux  and  exchange  terms  in  the 
energy  equation,  e.g.,  input  of  energy  by  wind,  contribution  from  mean 
flow,  work  done  against  pressure  forces,  etc.  This  would  provide, 
among  other  things,  a way  to  assess  the  effects  of  the  mean  zonal 
current  on  the  dynamics  of  the  system.  An  example  of  where  this  might 
be  particularly  useful  is  when  parameters  are  adjusted  so  that  an 
unstable  wave  will  develop.  In  this  case,  the  contribution  of  energy 
from  the  mean  flow  should  be  positive  throughout  the  integration. 
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A series  of  experiments,  similar  to  those  already  done,  but  using 
instead  a zonal  forcing,  should  result  in  odd-order  modes,  or  modes 
involving  anti-symmetric  v-components , including  Kelvin  waves. 

There  are  many  permutations  that  could  be  tried,  but  since  the 
model  is  highly  idealized,  emphasis  should  be  on  gross  aspects  of 
relevant  cases  and  not  on  differences  of  results  for  many  slightly 
different  sets  of  parameters.  It  is  felt  that  a problem  including 
continuous  vertical  dependence  should  be  studied  to  better  understand 
the  effects  of  atmospheric  forcing  in  the  light  of  vertical  energy 
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5.  CONCLUSIONS 

In  Chapter  2,  it  was  shown  that  there  is  evidence  for  the  exis- 
tence of  energetic,  several-week  oscillations  in  the  Atlantic  equatorial 
current  system.  While  it  is  difficult  to  determine  accurately  periods 
or  wavelengths,  it  appears  that  two  to  three  weeks  is  a reasonable 
estimate  of  the  former.  In  the  absence  of  any  definitive  observations, 
a longitudinal  wavelength  of  several  thousand  kilometers  is  indicated 
through  a rather  indirect  argument. 

Subsequent  analysis  of  current  velocity  data  reveals  that  advec- 
tive  effects  were  important  and  all  fields  were  highly  time-dependent 
during  GATE.  This  result  suggests  that  the  concept  of  a steady  equa- 
torial current  system  with  a corresponding  steady  meridional  flow  with 
upwelling  at  the  equator  may  have  to  be  modified,  and  that  it  may  be 
quite  difficult  to  determine  such  quantities  as  mean,  vertical  velocity 
observationally . 

The  presence  of  2-  to  3-week  wind  fluctuations  over  the  equatorial 
Atlantic  during  GATE  suggests  that  the  current  fluctuations  may  be  wind- 
excited.  Investigation  of  a simple,  two-layer,  baroclinic  system  shows 
that  oscillations  similar  to  those  observed  can  be  excited  by  impulsive 
meridional  forcing.  The  model  response  is  a resonance  phenomenon  rather 
than  a purely  passive  response  to  local  forcing.  That  is,  the  charac- 
teristics of  the  response  are  determined  chiefly  by  system  parameters, 
as  opposed  to  forcing  characteristics.  The  response  is  thus  composed 
of  the  free  modes  of  the  system.  In  the  absence  of  mean  currents,  the 
free  modes  excited  are  completely  consistent  with  the  theory  of 
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equatorially  trapped  waves,  and  furthermore,  most  of  the  energy  of  the 
response  is  found  in  the  several  gravest  modes  of  the  system.  The 
presence  of  mean  zonal  flows  tends  to  enhance  and  modify  the  response, 

I but  does  not  change  its  essentially  wavelike  character.  The  chief 

component  of  the  response  (termed  a Yanai  wave  in  the  case  without  mean 
flows')  appears  to  play  a role  similar  to  that  of  inertial  oscillations 
in  time-dependent,  mid-latitude  response  problems. 

While  the  presence  of  a free-mode  oscillation  in  the  equatorial 
ocean  could  be  the  result  of  a more  Indirect  forcing  (e.g,#  a wave  to 
wave  interaction  at  a meridional  boundary"),  wind-forcing  is  as  likely 
to  be  responsible  as  anything  else.  It  should  be  apparent  that,  even 
with  the  very  simple  model  described  in  Chapters  3 and  4,  the  combined 
effects  of  several  impulses  at  different  times  and  with  different 
wavelengths  would  result  in  a very  complex  mixture  of  waves.  It  is, 
therefore,  not  surprising  that  the  picture  emerging  from  GATE  is 
rather  complex. 
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